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ABSTRACT 


The expected motion characteristics of the real earth are sys- 
tematically analysed based on available dynamical theories for the rigid 
model, the elastic earth model and the earth model with liquid core. 
The various axes which are implicit in the dynamical theories are 
investigated regarding observability on the basis of astronomical obser- 
vations and suitability for defining reference directions. The observa- 
tional insignificance of the "diurnal polar motion" is demonstrated. A 
special effort ls made to clarify customarily used terminology. 
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LIST OF SYMBOLS 


This list of symbols gives definitions for the symbols used 
throughout the text. 


(X) E = (Xie, X 2E , X 3e ) t Geocentric and inertially oriented coordinate system. 

The third axis coincides with the north pole of the 
fixed ecliptic at a standard epoch. The first axis is 
fixed arbitrarily to the fixed ecliptic. 


(U) - (Ui, U2, U 3 ) t Body- fixed geocentric coordinate system. Its orien- 

tation within the body is given in a '’prescribed" way. 

(U)f The axes of this coordinate system coincide with the 

principal moments of inertia axes of the earth. 


(X)f 


(X)h 




• L 


H 


O) 


Same system as (U)f, but the first axis coincides 
with the node line of the fixed ecliptic and the 
equator of figure. 


The third axis coincides with the angular momentum 
axis (H). The first axis coincides with the node line 
of the fixed ecliptic and the plane orthogonal to (H) . 

Vector of luni- solar torque on the earth. 

Angular momentum vector. 

Instantaneous angular velocity of the earth. 


O 


Mean angular velocity of the earth. 


(I) 


Instantaneous rotation axis . 


(F) Axis of figure . 

(II) Axis of angular momentum. 
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(Eo) 

<C') 

(C") 

(C) 

L^. = (LiXp Ls Xp LaXp) 

Lx F ~ Li xp + i Eg x F 
(Uif, u 2F , 1 ) 

Up - UlF +iU2F 

Iu 

Ai , A 2 , A3 

CiJ 

U 

GMST 

k 

k s 


Axis of Eulerian pole of rotation. 

Axis of celestial pole for rigid and elastic earth 
model . 

Axis of celestial pole for the liquid core model. 

Axis of celestial pole . 

The components of L in the (X)f system. The same 
notation is used regarding other axes and coordinate 
systems. If no subscript is given, the components 
refer to the ( 0 ) -system. 

Complex form for equatorial component of torque L 
in (X)f system. The same notation is used regarding 
other axes and coordinate systems. If no subscripts 
are given, the components refer to hie (U)-system. 

Direction cosines of the axis (f) in the (U)-system 
whereby second -order terms in uif are neglected. 

Complex notation for the equatorial component of a 
unit vector along (F) in the (U)-system. 

Inertia tensor in (13)- system. 

Least, intermediary, and maximum moments of 
inertia. These are the elements in the diagonal of 
the inertia tensor if the coordinate system is selected 
such that the off-diagonal terms of the inertia tensor 
are zero. 

Small term in inertia tensor at position (i, 3 ) , 
Tide-generating potential. 

Greenwich mean sidereal time. 

Tidal effective Love number. 

Secular Love number . 
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(0 F , 4>f, <p F ) Euler angles relating the fixed ecliptic and the (U>- 

system. If no subscripts are given, the Euler angles 
refer to the (U)-system. tp is counted positive 
eastward. 

Xp Position of (F) in (X) E in complex notation. Similar 

notations are used for other axes . 

*Fh Free motion of (F) in (X) E . 

Xp f Forced motion of (F) in (X) E . 

x E p Precession of (F). 

Xp N Nutation of (F). 


Upf - UlFf + iu 1Ef 


UFh 

0pfj 0pf 
0FP, $FP 
0FN> $FN 

CTO 
O’ r 


Complex notation for equatorial components of unit 
vector of forced motion of (F) in (U)-system. 

Same as above except for free motion of (F). 

Forced motion of (F ) in obliquity and longitude . 

Precession of (F) in obliquity and longitude. 

Nutation of (F ) in obliquity and longitude . 

Frequency of Euler motion (rigid body). 

Frequency of Chandler motion (elastic body) . 



1. INTRODUCTION 


The strong increase in measurement accuracy, which is expected from 
the new laser generation, makes it necessary to look anew at the underlying 
principle of astronomical frames and their suitability for lunar laser ranging. 

No attempt is made in this study to introduce a new dynamical theory; rather, 
extensive use is made of available theories. Special efforts are made to 
identify the various axes, to discuss their relative body- fixed (with respect to 
the body) and space-fixed (with respect to inertial space) motions, and to 
investigate their dependence on defining constants, time varying parameters, 
suitability for providing a defining direction of reference, etc. 

The simplest theory is that of the rigid body earth model. The earliest 
investigations on this subject are several centuries old. The theory was 
essentially completed by Oppolzer [1882] when he included small second-order 
terms in his solution and thus derived the expressions for the diurnal body-fixed 
motions of the instantaneous rotation axis (I), which are a result of luni -solar 
attraction. Unfortunately, he re-emphasized an older concept, saying that the 
astronomical observations refer to that equatorial plane which is orthogonal 
to the instantaneous rotation axis (I) and neglected to investigate the question of 
observability anew in view of his expanded and more complete rigid earth theory. 
See [Oppolzer, 1882, p. 155]. Subsequently, the officially adopted set of nuta- 
tions has always been given for the instantaneous rotation axis (I) . It was only 
very recently that the question of observability was taken up again by Atkinson 
[1973, 1975, 1976] andOoe and Sasao [1974]. Severe difficulties of an habitual 
as well £is a practical nature (since a large amount of astronomical data has 
been reduced in a certain manner) surfaced at two recent international meetings 
at which this subject was discussed. In an attempt to correct the situation the 
General Assembly of the International Astronomical Union at Grenoble in 1976 



[IAU Transactions, 1977] passed a resolution to adopt the forced rigid body- 
nutations of the axis of figure (F) instead of those of the instantaneous rotation 
axis (I). However, in a subsequent meeting [IAU Symposium No. 78, 1977], at 
which a set for the nutations of the nonrigid earth was adopted, the resolution of 
Grenoble was ignored. It is thus appropriate to discuss the question of observa- 
bility in this study and to elaborate on the observational significance of the 
periodic diurnal body -fixed motions of the instantaneous rotation axis (I). The 
most accurate rigid earth theory available at this time is that by Woolard [1952] 
which was officially adopted. 

Besides the rigid model the motion characteristics of more elaborate 
models such as the elastic model and models with liquid core will be inves- 
tigated. The long history of research in the various earth models has not only 
produced more realistic models but also an amazing amount of confusion in 
terminology. In Table '1.1a summary of those phenomena which will be dis- 
cussed in great detail in this report is given including a summary of terms 
which can be found for them in the literature. Note that the terms in each box 
describe one and the same motion. No guarantee is given that the lists are 
complete. Table 1.1 is divided into two groups of motions. The body- fixed 
and space-fixed motions are called polar motion and nutation, respectively. 

This basic classification is the same as that given in [Munk and MacDonald, 

I960], The characteristic is that polar motion changes latitudes, whereas nuta- 
tion changes declination. The subdivisions are made according to the earth models 
and the cause of the motion. The force-free motions are similar in character to 
constants of integration and can exist independently of the forced motions It is 
also seen that the term "wobble" is associated with force-free body- fixed motions. 
The word "free" generally denotes force-free motions. In some cases even the 
body-fixed motions are called a "nutation. " Such terminology is used to illus- 
trate that any periodic motion should be called a nutation. If such conventions 
were tofbe followed, a unique identification of the motion components would 
require information as v to whether the motions were body- or ! space-fixed gnrl as 
to their periods . The'liquid' core model introduces new motrons in addition to 

if 
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Table 1 1 


Terminology Related to the Motions of the Instantaneous Rotation Axis (I) 


1# 

Model 

Excitation 

Polar Motion 
(Body- Fixed Motion) 

Nutation 

(Space-Fixed Motion) 

Rigid or Elastic 

1 force -free 

la. Chandler motion (elastic model) 

lb . Euler motion (rigid model) 

2. wobble 

free nutation 

luni- solar 
attraction 

1. Oppolzer terms 

2. forced dirunal motions 

3 dynamical variation of latitude 
(longitude) 

4. diurnal polar motion 

7 

1. astronomical nutations 

2. forced nutations 

3. nutations 

Liquid Core 

force -free 

1: nearly diurnal free wobble (NDFW) 

2 . nearly diurnal free polar motion 

3. nearly diurnal free nutation 

1 . assoc iated f re e nutation 
to NDFW 

2. free principal core 
nutation 

core resonance 

causes changes in diurnal polar motion and nutation (no name given) 




those of the rigid-body model. These are motions of the shell which result from 
interactions between the shell and the liquid core. The motions depend on the 
assumed structure of the shell and the core. If the core is assumed to consist 
of several layers of equal density then there exists the possibility of several 
nearly diurnal free wobbles. The free nutation due to the presence of the core 
is sometimes referred to as "free principle core nutation. " The word "core" 
should be interpreted in such a way that the core is responsible for the particular 
nutations of the shell . This motion, therefore, should be detectable from 
instruments located on the surface of the shell. Despite the many names in 
Table 1.1, the motions in each row are related by very simple kinematieal 
relations. According to the classical theory of Poinsot any rotational motion 
of a body around a fixed point can be represented by the rolling of a body cone 
onto a space cone. The line of contact is the instantaneous rotation axis (I). 

In Section 2.3.5, Poinsot's kinematieal representation is discussed in detail. 

It is not only the motion of the instantaneous rotation axis (I) which is of 
concern here. The following axes also have important body -fixed and space-fixed 
motion characteristics which will be investigated: 

(a) angular momentum axis (H) 

(b) axis of figure (P) 

(c) axis of Eulerian pole of rotation (Eo) 

(d) axis of celestial pole (C) 

The angular momentum axis (H) has the useful properly that its spatial motion is 
nearly independent of mass redistributions of the magnitude which is expected 
for realistic earth models [Eedorov, 1958]. The axis of figure (F) coincides 
at any instant with the direction of the maximum moment of inertia. The term 
"Eulerian pole of rotation" was introduced by Woolard [1952] . It corresponds to 
the position which the instantaneous rotation axis (I) would occupy if there were 
no external forces acting on the body. The term "Celestial Pole C" is used in 
this study to denote that pole to which astronomical observations and also lunar 
laser range observations refer. Its defining property is that it has neither 
periodic diurnal body -fixed nor space-fixed motions. It is now evident that the 
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expression "polar motion" is not unique; rather it is necessary to> name the axis 
to which it refers. 

The various motion components are explained most easily in terms of 
the corresponding mathematical expression. Therefore, a rather complete 
mathematical derivation is presented for the rigid and elastic models . Both 
models are treated together. With some simple specifications one can modify 
the results for the elastic body so as to obtain those of the rigid body. The 
derivations make heavy use of the works by Fedorov [1958] and McClure [1973] 
who gave a comprehensive derivation of diurnal polar motion. First, the nuta- 
tions of the axis of angular momentum (H) are computed. Since the spatial 
motions of (H) are nearly independent of mass redistributions inside the earth, 
the simple rigid earth model suffices for the computations. In fact, the deri- 
vations were done by Woolard [1952]. However, the most important steps in the 
derivation will be repeated using Doodson's [1921] tidal development in order to 
compute the limi-solar torques. This procedure can be found also in Melchior 
and Georis [1968] It has the advantage that it demonstrates the relationship 
between the nutations and the earth tides. In a second step the diurnal polar 
motions of the axes (H), (I) and (F) are derived and their observational signifi- 
cance is analyzed. The spatial motions of the axes (I), (F) and (C) will be 
obtained by transforming their diurnal polar motions relative to (H) into correc- 
tions which are to be added to the space-fixed motions of the angular momentum 
axis (H). 

The liquid core model is treated in a merely descriptive manner. Any 
rigorous mathematical treatment of such models is beyond the scope of this 
study. The various new motions are analyzed for their observability and the way 
in which they influence the choice of a reference axis. Probably the most classi- 
cal mathematical development for liquid core models in recent times is that of 
Jeffreys and Vicente [1957a, b]. 

The analysis regarding observability is based purely on geometric con- 
siderations. The observations are assumed to be reduced correctly (tidal 
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corrections, aberrations, etc.). In many cases the question of observability 
is analyzed conceptually, regardless of whether or not a particular measure- 
ment system is capable of reaching the required accuracy. In this report, 
the classical astronomical observations are analyzed. In [Leick, 1978] the 
analysis is extended to lunar laser ranging. ' 
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2. MOTION SPECTRUM POR RIGID AND 
ELASTIC EARTH MODEL 

2.1 Fundamental Differential Equations of Motion 


2.1.1 Axis of Rotation (I) 

Rl a geocentric inertial frame (X)e the motion relative to the center of 

mass is at every instant a rotation around an axis through the center of mass, 

in which the rate of change of the angular momentum vector H about the center 

— > 

of mass is equal to the resultant torque L. 



The time derivatives refer to components on the inertial axes since this form of 
the equation of motion holds only in an inertial system. 

The equation of motion involving the time derivatives of components 
with respect to the moving axes (U) is [Goldstein, 1965] 




+ O) X H 


( 2 . 1 - 2 ) 


CO is the angular velocity vector of the moving frame with respect to the inertial 
frame . The angular velocity vector lies along the axis of infinitesimal rotation , 
a direction which is also called the instantaneous axis of rotation (I) . 

The equations (2.1-2) are a system of three linear first-order differen- 
tial equations They are usually referred to as Euler's dynamic equations . 
Substituting equation (2. 1-1) in equation (2.1-2) and making use of the usual 
summation convention, which calls for a summation over repeated indices, 
equation (2. 1-2) is equivalently written as 


Li 


dHi 

dt 


£ ijk ^ 5 Hk 


ij j > k 1,2,3 


(2.1-3) 
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All components in the equation above refer to the moving axes . e ijk is the 
"alternating" tensor, defined by the following properties: 

0 if any two subscripts are equal 
Cijic = 1 for even permutation 

-1 for odd permutation 

Equations (2.1-3) are quite general. They represent the rotational motion for a 
rigid or elastic body. The standard expressions for the angular momentum, the 
tensor of inertia, and the relative angular momentum are, respectively 

Hi = Iij Wj + hi 

lii = / (U k Uic 6ij - U* Uj) dm 

m (2 . 1-4) 

hi = f ^ijk Uj Uk dm 

The symbol 6ij above denotes the Kronecker delta. The relative angular mo- 
mentum is caused by internal motions of the mass particles themselves. The 
dot in Uk denotes differentiation of Uk with respect to time. The integration is 
taken over the whole body of mass M. 

Substituting equations (2. 1-4) in (2.1-3) leads to yet another form of the 
equations of motion, usually referred to as the Liouville equation 

Li = ( Itj w j + hi) + f ijk aJ j ( Li + hk) (2.1-5) 

The orientation of the rotating coordinate system (U) is arbitrary at this point. 
An extensive discussion on the possible choices can be found in [Munk and 
MacDonald, I960]. Obvious candidate axes are the principal axes for which 
the products of inertia are zero. The body-fixed system (U) to be used in this 
study 'is attached in a "prescribed way" to some observatories, and its origin 
is at the center of mass. Thus the (U) -system and the axis of figure system 
(U) F do not necessarily coincide. The introduction of the (U)-system allows a 
slightly more general formulation of polar motion inasmuch as it introduces 
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constant components in tlie off-diagonal positions in the moment of inertia 
tensor . These terms are responsible for a constant polar motion component 

t 

which is equal to the angle between the third axis of the (U)-system and the 
axis of maximum moment of inertia. The periodic polar motion components 
are independent of the specific choice of the (U) -system This fact justifies the 
use of the term inertia "tensor." It is thus permissible to select the coordinate 
system according to its advantages in formulating a specific motion most 
simply. At this stage of the development it is only required that the third axis 
of the system (U) be nearly aligned with the rotation axis such that the products 
and squares of certain small quantities are negligible . 

For evaluating the Liouville equation, the inertia tensor has to be known. 
In order to emphasize that the constant components are included in the off- 
diagonal terms, the subscript U is used. Thus, 


Iu 


'Ai + cu 
Cl2 
Cl3 


C33 

Ai+ Css 
C 33 


Cl3 

023 

A3+ C33 


(2.1-6) 


Ai and A 3 are the least and the maximum moments of inertia. The Cij are small 
quantities compared to Ai and A3. The angular velocity components in (U) are 



■wr 


Uii ‘ 

00u - 

C0 3 

= 0 

U21 


_ w 3 _ 


_1+U31_ 


Q, is the mean angular velocity of the earth rotation The dimensionless numbers 
u« are very smalL. In fact, uii as well as cij/As and hs/OA 3 which appear in the 
Liouville equation are of the order 10 6 or even smaller for realistic earth 
models . Neglecting second-order terms in these quantities gives a first-order 
theory accurate to order 10~ 6 . Equation (2.1-5) reduces to 

Li — Ai Uii^ + Ci 3 £2 + hi + U 21 (A3 - Ai) - C23^1 - h 2 ^ 

L 2 = Ai U 21 £2 + C23^1 + I 12 - Uii (A 3 -Ai) Cl 2 + ci 3 £l 3 + hi£l (2.1-8) 

L3 — C33 + A3 U31 C 2 + I13 
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u XI and us! appear only in the first two equations of (2. 1-8), whereas U31 is only 
present in the third equation. This makes it possible to separately solve for 
the quantities (uu, u 2i ) and (U31.), respectively. Using complex notation, the 
first- two equations in (2 1-8)' are combined as 

u r = iffr(ui-^) ( 2 . 1 - 9 ) 


The symbols have the following meanings: 

i = J~L 

ui = Uu + iuai 

cx r = a 

Ai 

^ iL | c ic t h ih 

(A3 -Ai)fi 2 A3-A1 (A3 -Ai)^ (As-Ai)^ (A3 - Ai)£^ a 


The complex quantities are: 


+ i^2 
L = Li + i L2 
e = Cj3 + i C23 
h = hi + ihs 


The dimensionless function ^ is called the excitation function. It is a function 
of the torque, the inertia tensor, and the relative angular momentum . The 
third equation in (2.1-8) becomes 


• _ L3 _ C33 _ I13 

A3 £1 A3 A3 ^ 


( 2 . 1 - 10 ) 


The complex variable ux completely describes the body-fixed motion of the 
instantaneous rotation axis (I) in the -frame (U). The direction cosines of (I) 
are obtained from equation (2.1-7) as 

= (ail, usx, 1) (2.1-11) 

M 

where second-order small quantities are omitted. 
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« i 2.1.2 Axis of Angular Momentum (H) 

>'The motion of the angular-momentum vector in the (U)- system is given 
■ by equation (2. 1-4). Neglecting,’ once again, the second-order terms, equation 
(2. 1-4) becomes 1 • '• 

H = Aifi uj + O c + h 

* H3 “ A 3 ^i (l+Ual) + C33O + ll 3 ^ 

with 

H = Hi + iH 2 


Thus, substituting the solution of equation (2. 1-9) for the instantaneous rotation 
axis (I) into equation (2. 1-12) gives the body-fixed motion of the angular momen- 
tum in the (U)-system. The direction cosines, accurate to the first order, are 

*) s < Ul »' U3 ” x > < 2 - x - 13 > 

t 

The position of the angular momentum vector in the terrestrial system is 
largely a function of the position of the instantaneous rotation axis which is given 
by ur. But the relative position of these two axes is, in addition, directly 
dependent on the disturbances of the tensor of inertia' and the relative angular 
momentum. Equations (2.1-12) and (2.1-13) give 



u H - Uj 


A3 -Ai 
Ax 



h 

A3 £2 


An estimate for the coefficient of uj on the right-hand side is 


A3 - Ai 

Ai 


0.0033 


For an elastic earth model or in any other nonrigid model for which c and h 
are not zero, a more significant separation is expected between the angular 
momentum ‘axis (HJ and the instantaneous rotation axis (I). 

, . ■ . ~ t 
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2.1.3 Axis of Figure (F) 

The motion of the axis of figure is a function of the disturbances of the 
inertia tensor. In the principal axis system (U)f, in which the third axis coin- 
cides with the axis of figure (maximum moment of inertia axis) , the off-diagonal 
elements of the inertia tensor are zero. The third axis of the two systems (U> 
and (U) are related by two small orthogonal rotations as follows: 

(U)f = R 2 (uif) Ri (-u 2 f) (U) = R (Utf, U2f) (U) 


Uif and u 3 f are the direction cosines of the axis of figure (F) in the terrestrial 
system (U). Next, the inertia tensor I u of equation (2.1-6) is transformed to the 
principal axis system. The transformation properties of the components of the 
inertia tensor are determined by the fact that the matrix Iu transforms under R 
by a similarity transformation [Goldstein, 1965, p. 147] . 

If = R I u R 1 


To the first order, the inertia tensor in the (U)f system is, therefore, 


If - 


Ai + cii 

Cl2 

[Uif(Ai -A3) + C 13 ] 


Cl2 

Ai + C 22 

[U2F(Ai -A3) + C23] 


[ujlf(Ai -A 3 ) + C13] 
[Usf(Ai-A 3 ) + C23] 
A3 + C33 


(2.1-14) 


The diagonal elements remain unchanged after the transformation. Since the 
off-diagonal elements are zero in If, at all times, the following relations for 
the direction cosines are taken from (2 . 1-14): 


Uif 


Cl3 

A 3 _ Ai 


U 3 f 


CS3 

A3 — Ai 


(2.1-15) 


The small terms cn, cia and C23 which are present in If do not enter the 
excitation function S? explicitly. - In-case.cis 5 ^ 0, an additional rotation around 
the third axis will complete the diagonalization of the inertia tensor. Such a 
rotation, however, is not of concern here when treating only the aspect of polar 
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motion. Equations (2.1-15) clearly show how a constant term in C 13 and C 23 
causes a constant component in Uif, i. e. , a constant offset between the third 
axis in the (U)-system and the axis of figure. In complex notation the direction 
cosines are written as 

Up = UlF + iu 2 F 

With the help of (2 . 1-9), equations (2. 1-15) which determine the direction of the 
axis of figure (F) in the terrestrial system (U) are written as 

„ _ o (2.1-16) 


2.1 4 Coplanar Motion 

Combining the equations (2.1-11), (2, 1-12), (2.1-13) and (2. 1-16) 
gives an interesting relation: 

u H - up = + ^ (2.1-17) 

The presence of the relative angular momentum in equation (2. 1-17) is some- 
what disturbing in view of the geometrical interpretation. If 

h = 0 (2.1-18) 


then the angular momentum vector, the axis of figure, and the instantaneous 
rotation axis are located in one plane, as is expressed by 

Ai 

UH - Up = — (U! - Up) 



Such a motion is certainly true for a rigid body motion where h = c = 0. The 
separation between the various axes is a function of the moments of inertia. 
For a spherical body with Ai = A3, the separation between (I) and (H) vanishes . 
Since for a realistic earth model the ratio A1/A3 is smaller than unity, the 
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angular momentum axis lies always between the axis of figure and the instan- 
taneous axis of rotation. These motion characteristics are valid whether or 
not external torques act on the b.ody. 

2.2 Spatial Motions of the Axis of Angular Momentum 

2.2.1 Euler’s Kinematic Equations 

The relationship between the coordinate system (U) and the inertial 
system (X) E , i. e. , the ecliptic system at a standard epoch To, is given by 
Euler’s kinematic equations. The situation is demonstrated in Figure 2.1. 



Figure 2. 1 Euler’s Kinematical Relations 
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The Euler angles (6, 0 , 9) and their time derivatives ( 0 , 0 , 9) repre- 
sent the motion of the (U)-s’ystem in space . But since this system is rigidly 
connected to the body, the angular velocity components ah as per equation 
(2. 1-7) describe the same motion. The relation between the two sets of angular 
velocities is found by resolving each of the angular velocities 0, 0, 9 along the 
{U l5 U 3 , U 3 )-axes and adding the components along each axis as follows: 

<*>! = - 0 cos 9 - 0 sin 0 sin<p 

co 2 = 0sin9 - 0 sin 0 cos 9 (2.2-1) 

co 3 = 0 cos 0+9 

The inverse relations are 

0 sin 0 = -t 0 isin<p -co 2 cos 9 

0 - - co i cos 9 + co s sin 9 (2.2-2) 

9 = co 3 - 0 cos 0 

2.2.2 Poisson's Equation of Motion 

The spatial motion of the angular momentum vector is the least sensitive 
axis regarding internal mass movements. Mass redistributions affect the 
motion of (H) in space only through their effect on the luni-solar torques. 
Fedorov [1958] made an extensive analysis of these disturbances. His calcula- 
tion showed that the luni-solar torques as computed for a rigid and perfectly 
elastic earth model differ only in the order of 10 -s . Within such an accuracy 
level, therefore, the computation of the nutation of the angular momentum can 
be based on a rigid earth model. 

Here a rigid model with equal least and intermediary moment of inertia 
is selected. Although the system (U) was used in Figure 2. 1 for the purpose of 
a more general representation, permitting a separation between the Ua-axis 
and the axis of figure, the following derivation refers to the (X) F -system, whose 
third axis coincides with the axis of figure (F) and whose first axis is along the 
node line as defined by the fixed ecliptic and the equator of figure. There should 
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be no confusion about the various coordinate systems used here. The (X> or 
(U> -systems are especially useful for dynamical purposes since the tensor of 
inertia is diagonal -in this, system. (U)f and (U) are related via equation (2.1-16). 
A further simplification is achieved by assuming the same magnitude for the 
least and intermediary moment of inertia, i. e. , Ai = As. In such a case the 
selection of the first axis in (U)f is arbitrary. 

The coordinate system (X) F does not take part in the daily rotation of the 
earth. Its orientation in space changes only due to precession and nutation, 
which is being expressed by the small rotational velocity oo . The components of 
oo 'onthe axes of (X) F ar e obtainedby substituting <p~0 and 0=0 in equations (2 . 2-1) : 

CO\ = - Qf 

00 2 ' = _ j^> F sin0F (2.2-3) 

co 3 ' — $f cos 6p 


These velocities can be substituted in Euler's dynamic equations (2. 1-3). 
The angular momentum is 


Hi = Ai to t i = 1 , 2 

H 3 - As(C^3 + cp) = A 3 £03 


(2.2-4) 


The velocities 00 / in these equations are due to the motion of the frame (X) F with 
respect to (X)e, and co 3 is the velocity component of the earth rotation about the 
third axis with respect to the same system . Equations (2.2-4) are valid at any 
instant. Since the moments of inertia about all the axes perpendicular to the 
axis of symmetry are the same for the model to be considered here,, the moments 
of inertia Ai do not change with time even though they are not referred to an axis 
attached to the body. 

Euler's dynamic equations (2.1-3) in the (X) F system are 

A1OC1 + 0O3 A3 (£03^+ <P) - co 3 / Ai£^2 — Lix f 

A ]£h 2 - £^iA3(£d3 + < P) - £0 3 / Ai£0 i = Latp (2.2-5) 

A3(£^3+<P) + £ 0 iAi£ 03 - CO 2 A 1 OO 1 = L3x f 
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These equations are rigorous . Their solution, in combination with (2. 2-3) 
would yield the nutation of the axis of figure. No attempt is made to solve this 
system of equations, rather the motions are investigated which result if 
certain simplifications are introduced, such as neglecting the product of small 
quantities Oil and their derivatives. The first two equations of (2 ..2-5) then 
give 


co P = 


LiXp 

A 3 ( P 
Tex, 


( 2 . 2 - 6 ) 


Oil = - 


Aa<P 


The third component Iex F is identical to zero because of symmetry relative to 

> 

the equator (compare Section 2. 2.3). Therefore, thethirdequationin(2. 2-5) gives 

AaP = 0 

which has the solution 

ip = constant 

ip is the mean velocity of the earth rotation. Previously we denoted this mean 
velocity by £2 so that we have the identity in notation 

<P = Q, (2.2-7) 


Using equations (2.2-3) and (2.2-7), we obtain from (2.2-6) the well-known 
Poisson equations of motion 


0 sin 6 


0 


Lixp 

“ Aafi 
L2 X f 
A3 £2 


( 2 . 2 - 8 ) 


These equations are of fundamental importance for describing the orientation 
of the earth. Indeed, the Poisson equations are the virtually rigorous equations 
for the nutation of the angular momentum axis (H). To show this, consider the 
coordinate system (X)H - in which the third axis coincides with (H) and the first 
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axis is along the node line of the fixed eliptic and the plane orthogonal to (H). 
Introduce the Euler angles $h and 0« in order to relate the two systems by 

(X) E = Rs(-0h) Ri(9h) (X)h (2.2-9) 

The components of H in the (X) E system are 

- sin$H sin 0 h 
cos 0 h sin 0 h 
cos 0 h 

where | H| = A 3 £1. Equation (2.1-1) can now be solved directly for the motion of 
the angular momentum axis in space. The first two equations are 

” ("A^ sin$H sin0 H ) = Lix e 

^ (A 3 ^ cos^h sin0H) = L 2 x e 

Lix E an( i Lsx E are torque components in the ecliptic system (X) E . The two 
equations above are solved for the rates of the Euler angles . 

0 h = (L 2 x e cos$h - Lix E sin$ H )/(A 3 fi cos0 H ) 

( 2 . 2 - 10 ) 

(£> h = (-Lix e c °s^h - L 2 x e sin^H)/(A 3 ^ sin0H) 

The torque in the (X) h system is related to the torque in (X) E according to 
(2. 2-9) as 

L Xh = Ri (-04) Rb (0h) Lx e (2.2-11) 

The difference in L Xu and L X|r is negligible for all practical purposes because the 

n * 

systems (X) F and (X)h are closely aligned. In anticipation of later discussions on 
the free and forced motions , it is noted that in the present context, only the 
effect due to the difference in forced position enters. With this approximation 
and knowing that L 3 x F = 0) equation (2.2-11) gives 

Lix f = Iax E cos$h + L 2x e sinJ/>H 

L 2 x f . = (-Lix E sin^) H + Lsx e cos!/>h)/cos0h 
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Solving these equations for Lix e and L 3 x E , and substituting them in equations 
(2. 2-10) results in the Poisson equations (2. 2-8). 

The most thorough solution of Poisson’s equations can be found in 
[-Woolard, 1952]. He made use of a second-order force function for the gravita- 
tional attraction. The positions of the moon and the sun were taken from the 
respective theories of Brown and Newcomb. Woolard’s form of the Poisson 
equation is his equation (30). The verification of identity between his form and 
equations (2.2-10) can be made almost by inspection. To be more' explicit, take 
the right-hand side of Woolard’s equation (4) which gives the torque components 
in the (U)f -system as follows: 


* 


r 3 u 

Li 


ae 

Lz 

= - R 3 (<P) 

af /sm6 

Ii 3 

Uf 

0 


where U is his force function. The third torque component is again zero 
because of symmetry . 

The torque components in the (X)f system are 


or 


Lx f - Ha (-<P) 


Lx 

L s 

Ij 3 


J Uf 



-N ^ 

3U 

“ ae 



/ sin 0 


0 


Substituting the equation above in (2.2-8) gives exactly Woolard's form of the 
Poisson equations . 


19 



Historically, the Poisson equations have played an important role in 
deriving a nutation set. The presently adopted set of nutations is that of the 
angular momentum axis (H) based on Poisson's solution [AENA Supplement, 

196-1, p. 44]. Off ieiallythe adopted set is termed the ^nutation of the rotation 
axis" based on a rigid earth model. But as will be demonstrated in great detail 
in later sections of this study, the separation between the angular momentum 
axis (H) and the rotation axis (I) based on a rigid earth model is smaller than 
0'.'0002, which is the accuracy with which the adopted set of nutations are given. 
The fact that the angular momentum axis (Poisson solution) is a good approxima- 
tion to the rotation axis can readily be understood by considering equation (2.1-19). 
This relationship was first pointed out by Oppolzer [1882] . In fact, in earlier 
times one dealt only with Poisson's equation when computing the nutation and, 
therefore, one was not able to discover the nearly diurnal polar motion terms 
of the instantaneous rotation axis (I) . It was Oppolzer [1882] who pointed out 
the importance of the small terms which wereneglected whenderivingthePoisson 
equations from the dynamical equations (2.2-5). He included these small terms 
in hi's calculation and demonstrated that the rotation axis and the angular momen- 
tum must have a diurnal periodic body-fixed motion. The expressions which 
describe these periodic motions in a body-fixed frame are since called the 
"Oppolzer terms" [Woolard, 1952, p. 159], The corresponding expressions 
which have to be added to the nutations of the Poisson solution in order to arrive 
at the nutations of the axis of figure (F) or axis of rotation (I) are unfortunately 
called "diurnal nutations, " although they do not have any diurnal period whatso- 
ever [Woolard, 1952, p. 132]. In Section 2.3 a detailed representation is given 
for all phenomena mentioned above. 

2. 2. 3 Development of External Torques 

The gravitational potential of the earth at a distant point D (re , $ D , A 0 ) 
of unit mass is [Heiskanen and Moritz, 1966]: 


20 



V = 


GM 

r 0 


-sm' 


[<Jnm Ena (^0 >A.d) Knm Sum (^Dj -^- d)] ( (2.2~12) 


The symbols have the following meanings: 

# 0 , Ad latitude, longitude of the distant body in the axis of figure 
system (U)f 


ro 


geocentric distance ot the distant body 
mean earth radius 


G constant of gravitation 

M mass of the earth 

Jnn » Knm potential coefficients 


The functions K and S are 

Ad) = T? aa (cos ®d) COS m A D 
Sno (*Df Ad) = Pnm (cos $ D ) sin m A D 

s 

where P na are associated Legendre fu nctio ns . 

Since the origin ot the coordinate system (U)f is at the center of mass of 
the earth, the first-order potential coefficients are zero: 



Jio — Jii — Kao = 0 

Using this condition one can rewrite equation (2.2-12) as follows: 



J n P„ (sin <?□) 


(Sin ®o) Pnta t Knm Sum] 


The coefficient J s is approximately l(f times larger than any of the other coeffi- 
cients . The second-order approximation of the potential is 
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v ^ ^ r"fej j - p 3( sin *>>J 


This equation.can be modified by using the well-known relations 

Ai + Aa 


A 3 


J 2 


M a° 


(2.2-13) 


and 


as 


P=( Sto *„) = 


or 


V = f(r D ) 


V = f(r D ) 


3GM 2 T TT s 
~2? r a J2 U3F d 


3G 

2r 


6 (As - Al) U 3 f* 


(2.2-14) 


(2.2-15) 


It is assumed, as usual, that the least and intermediary moments of inertia are 
equal, i. e. , Ai = Aa. U3f d is the third coordinate of the distant body in the (U)f- 
system. The potential field f(r D ) is a central field; it has no effect on the torques. 
The expression of the potential (2.2-15) is identical to the one used by Woolard. 

The general expression for the torques is now readily obtainable. The 
gravitational force between the element dM and the mass of the disturbing body 
Mo, taken as a point mass, according to Figure 2.2, is 

dF = GM d dM -z - - M d grad (dV) 
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Figure 2. 2 Luni -Solar Torque 



Substituting this expression in (2.2-16) one obtains the torque components 
along the axes of (U) F as 
3M D GMa 2 J s 

Lu f ~ r 0 5 (Y d Zo , -XdZ 0 , O) (2.2-17) 
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The third torque component is zero. The total torque is orthogonal to the 
earth-moon line, and zero when the moon crosses the equator. Evaluating 
the torques requires knowledge of the time-dependent positions of the disturb- 
ing body in the 'system (U)f . Such positions are implicit in Doodson 1 s [1921] 
tidal development. Doodson's development is used here not only because it 
gives the positions of the disturbing body in the system (U) f , but mainly 
because it is anticipated to express the torques, and thus the nutations, in 
terms of harmonic series containing the tidal frequencies. This procedure 
will make it possible to study the relations between nutation and tidal theory. 

In order to make the comparison in the subsequent development 
easier the Cartesian coordinates in (2.2-17) are replaced by spherical 
coordinates. Using the well-known relation 

Psi (sin<Po) = 3 sin<pb cos<Pd 
the torque in the (U) f system is 

Q £ 'J 1 

Lu f = -*9 MM D a J 2 [sinAfj P 3 i (sin<p D ), - cosAo P 2 i(sincpo), O] 

r ° (2.2-18) 

Doodson’s [1921] expansion of the tide generating potential is now 
given in the notation of McClure [1973] : 

U X! (“ ) (sin<p)y'A nB j D cos [a^ + mA+ (n-m)-^] 

fr 1 c & n=2 m=0 \ Cd / f-' 

D (2.2-19) 

The tesseral coefficients A S1 jD are related to Doodson's original coefficients 

by _ _ 

A 2 ijd = - , S'A 2 ijD (2.2-20) 

The index j denotes individual tabular entries. The other symbols in 
(2.2-19) have the following meaning: 
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- D ■ • .disturbing-.body (sun, . moon) 

,Md • ’>mass of disturbing body 

Cd * . • - mean .distance- of disturbing body from, geo center 
G ~ , constant of gravitation. - 

a mean earth radius 

G a a i d tidal argument of body D 

(f> , A latitude and longitude of a station in the axis of 

>’ > -t 

figure system (U) F 
TT constant 3.14. . . 

The tidal argument can be computed from the relation [McClure, 1973] 

■ «nmj + mA= di T+ (cfe-5) s + ,(d3 - 5)h 

+ (d4-5)p+ (ds-5)N / + {d s -5)p 0 (2.2-21) 

b 

where all di’s are integers. For diurnal tide components di is equal to one. 
The other symbols on the right hand side of (2.2-21) are Doodson’s 
standard variables, which are sometimes referred to as the mean longitudes. 
They are related to Brown’s fundamental arguments ( 1 , 1* F,D, Q)as 

mean longitude of the moon 
mean longitude of the sun 
longitude of lunar perigee 
longitude of ascending 
node of the moon(£}) 
longitude of the perihelion 

The mean longitudes are measured from the mean vernal equinox. Expres- 
sions for Brown's fundamental arguments are given in AENA Supplement 
[1961, p. 44] in terms of polynomials of time. T is the local mean lunar 
time; it is reckoned from the lower transit of the moon. F inally , the Green- 
wich mean sidereal time is [Doodson, 1921] 

GMST = T + s- ir_A ’ (2.2-23) 


s = F + £2 
h = F-D + a 

p = 1 + F + £2 

( 2 . 2 - 22 ) 

N'=-£2 

P s “- l'+ F - D + £) 
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This relation assumes that the zero longitude is at Greenwich. Such an identi- 
fication is possible because the least and intermediary moments of inertia 
are equal in our basic dynamical model. For reasons of abbreviation the 
tidal argument (2.2-21) is sometimes given in the code form 

di ds cb. d^ ds d6 . 


As can be verified with the help of equations (2.2-21) and (2. 2-23), the code 
form, e.g., 165.555, denotes the tidal argument with sidereal frequency, 
which sometimes is also symbolized by S Ki or 111 Ki depending whether it 
results from solar or lunar attraction. The positions of the disturbing body 
in (U)f are obtained by comparing a formal spherical harmonic expansion 
of the tide generating potential with equation (2.2-19). Such a formal expan- 
sion is [Melchior, 1966, p. 15], 

U== E^“E t np (oosy) (2 - 2 ' 24) 

D n=2 

The angle y is the spherical distance between the station vector and the 
vector to the disturbing body, so that 


cos 7- sin sin <Po + cos <f> cos<Pd cos (A-A 0 ) 


Using the decomposition formula (Heiskanen and Moritz, 1966, p. 33) the tide 
generating potential (2.2-24) is written in the general form 


U 


co n 

V- GM D y V' 
^ Co n^S^O 


a_ 

r 0 


n 


W„ 


n m 


p. 


(sin<p) Pm, (sin<£b)cos m(A-A D ), '2.2-25) 


where 



2(n-m) ! 
(m-m) ! 


and 


Wcu = 


i 
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Substituting 

cos m(ArA 0 )= cos mAcos m A D + sin m A sin mA p 

into equation (2.2-25) and replacing the tidal argument in equation (2.2-19) by 

<< - - 
ft 

. cos (a n8 + mA+ (n-m) -) 


= cos(mA)cos a aBjD + (n-m) ^ -sm(mA)sina nBii0 


+ (n-m) J 


a comparison can be made of the coefficients of sinm Aand cos m A in the 
equations (2.2-25) and (2.2-19). The following relations are found: 
n+1 



V&m (sin<po) 


cos m A 0 
sin m A d 



cos 

-sin 


Kmjo + (n-m) f 
r c 

Uwjo + ( n ~m) v ~ 


(2.2-26) 


Only the terms with‘n = 2 and m = 1 enter into the torque equations 
(2.2-18). Combining equations (2.2-18) and (2.2-26) finally gives 
expressions for the torques 

Li Up = GMM d 3J 2 ^J A 21 j d cos (O SU0 ) 

j 

a 2 v— ' 

Lau F = GMM D ~ 3 J a / . As^p sin (a 2 ij 0 ) (2.2-27) 

c D j 

L3u f = 0 

} » , J, i 

If J 2 is replaced by (2.2-13) and the equation (2.2-20) is used then the 
equatorial components of the torque in the (U) F -system are in complex 
form: 

L Uf = 2 Aj 0 e ~ l(ajo) (2.2-28) 

j 


27 



where 


A J0 = f ^'(As-Ai) A 3Uo (2.2-29) 

The subscripts 2 and 1 are omitted in the exponent of (2.2-28). As an 
intermediary result it is seen that the torques and thus the nutations 
depend only on the tesseral harmonic coefficients A 31 . The torque com- 
ponents in the system (X) F are 


Lx r - L Uf e 


up 


A ^ e 
j 


-i(G!j D -<p) 


(2.2-30) 


c p is the Euler angle which measures the earth rotation. 

The expressions (2.2-30) can be modified to express the torque in 
terms of frequencies which are symmetric with respect to the earth rota- 
tion (sidereal frequency). Combining equations (2.2-21) and (2.2-23) and 
taking only the second order tesseral, i.e., n = 2 and m = d t = 1, gives 
the tidal argument 


- GMST + IT + (ds — 6)s + (d3-5)h + (dt-5)p 
+ (d6-5)N» + <de-5)p s 


(2.2-31) 


Combining the latter terms this tidal argument can be rewritten as 

a s = GMST + V + Affj (2.2-32) 

where A eg is called the nutation argument. The mean sidereal frequen- 
cy is 

Q= (GMST) ’ (2.2-33) 
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Two frequencies are symmetric to fl if 


Oi\ 4 . = £2 + Attj-j- 


ttj- = £2 + Aaj_ 

(2.2-34) 

and 


Aa,. = -A aj + 


hold. The positive sign in A dh + denotes that a* > £2. 
tidal arguments the symmetry is expressed as 

In terms of the 

a J + = GMST + it + Agj + 

CHj _ = GMST +’ 7T + A ctj _ 

(2.2-35) 


Such a splitting up of tidal frequencies is implicit in Doodson's develop- 
ment. The respective symmetric frequencies in the expression (2.2-30) 
can now be combined with the help of equations (2.2-35) as follows: 

Ai+e -m + -v) + Aj _ e -m*, --** 


. . D -i(GMST-<p+w-Ao' J+ ) A . -i(GMST-<pftr-Aa,+) 

— aj + e + Aj _e 

(2.2-36) 

= e ^ [ - (A, + + Aj _)cos Attj + + i (Aj + - Aj _) sin A ctj + ] 

The small difference between the Greenwich mean sidereal time, GMST, 
and the Euler angle is due to luni-solar nutation and planetary preces- 
sion. It causes the coefficients in the expansion of the torque to have a 
minor time dependent variation. The most important of these is due to 
planetary precession because it is secular. It is responsible for sm all 
secular terms which can be found in die nutation tables [Woolard, 1952, p. 
153] . For the purpose of discussions in this study, the approximation 

GMST = <p (2.2-37) 

is made whenever torques are computed. It may be emphasized that the addi- 
tional difference between GMST and <p which is due to the free motion is not 
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important in calculating the torques. In fact, all our derivations, so far 
as they are dealing with the forced motion, correctly assume that polar 
motion is zero. Polar motion, a result of the homogeneous solution, is 
independent of the forced solution. 

Prom equations (2.2-30) and (2.2-36) the expression for the torque 
becomes 

L Xp 1“ (A -3 + + A 0 ~) cos Aa j+ + i(Aj+ - A ,-) sma 3 + ] (2.2-38) 

3 

Thus, two tidal waves which are symmetric to the diurnal frequency form 
only one constituent in the torque. It is understood that the summation 
includes both the lunar and solar terms contained in Doodson's tidal 
expansion. 


2.2.4 Integration of Poisson's Equations 

The solution of the Poisson equations gives the nutation of the 
angular momentum axis (H). Denoting the derivatives in the Euler angles 

by 

♦ 

K = Q H + i^H sin 0 (2.2-39) 

then from the Poisson's equations (2.2-8) it follows 
i L Xp 
= " a As 

The homogeneous solution (force free, L = 0) is 

rv 

* Hl1 c • (2.2-40) 

where c is a complex constant. This expresses the simple fact that the 
angular momentum vector forms an invariant direction in space if no external 
forces act on the body. 
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/ 

A particular solution ( forced motion, L ^ 0) is given by 

** - ^r/ L(T)dT 


(2.2-41) 


= E~( A J+ +A J-) cos Aaj + +i(A j+ -Aj_)sinAa J+ ] dr 


where the luni-solar torque of equation (2.2-38) is substituted. A«j is f 'a 
polynomial in time. Before integrating, sidereal terms with = 0 are 
separated. The sidereal terms have per equation (2.2-31) the argument 
number 165.555. They cause the secular terms in the expressions of the 
Euler angles and constitute the luni-solar precession. There are two 
waves with sidereal frequency called a K : , and 8 Ki , which are due to 
lunar and solar attraction, respectively. Thus, e.g. (2.2-41) gives for 
precession and nutation 


i — ^Hp + Hh n 


T^T" (A* + A s )t 

A 3 U ki «i 


(2.2-42) 


+ 


1 


“(At+-At _) 


A 3 O Aa j + 


cos + i 


(A. J + + A 4 _) 


A ttj 


+ 


sin Acc j + | 


Using the torque in the form (2.2-30) and substituting equations (2.2-31) 
and (2.2-37) an equivalent expression for the nutations can be found 


X HN 


_L..y^-A.i.. e -iAftj 

Aaflii-'Affij 

j 


(2.2-43) 


The summation goes over all j , excluding those for which A aj = 0 . 
The luni-solar precession x HP of the angular momentum axis (H) is 
according to equations (2.2-39) and (2.2-42) 


0 HP = 


A 3 sin 9 q 


(Am K + Ab k ) 


0 H p = 00 


(2.2-44) 
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The two tidal amplitudes of equation (2.'2-44) are related to the actual 
values given in Doodson's table by equation (2.2-29). Doodson's 
coefficients are negative so that one obtains a circular westward motion 
of the angular- momentum axis due to precession. The subscript HP 
indicates that the elements of precession refer to the angular momentum 
(H) . This is especially important for the obliquity 9 H p • It will be 
demonstrated later that the obliquity of the celestial equator (Section 2.3.4) has 
a slightly different value. Since 6 Hp is a constant it appears as if there 
is no luni-solar precession in obliquity, i.e., the fixed ecliptic and the 
mean equator always subtend the same angle. But this is only true 
because the planetary precession was neglected when computing the torques, 
as is implied by the approximation (2.2-37). The inclusion of the planetary 
precession results in small non-linear precessional terms in both longi- 
tude and obliquity. A more complete expression for luni-solar preces- 
sion is 


0n P = fit + f 3 t 2 + 
0 Hp = So + 6 at 3 + * • • 


(2.2-45) 


fj, is called the constant of luni-solar precession and can be identified 
with the corresponding factor in equation (2.2-44). The coefficients fs and 
Q a are computable from theory. The respective expressions are given in 
[Woolard and Clemenee, 1966, p.242]. Note that there is no linear term, in ob- 
liquity. For the sake of completeness it is mentioned that the use of 0 O instead 
of 9 in the first equations of (2. 2-44) and (2.2-45) implies that the nutations Xhn 
should be amended by very small periodic terms whose arguments are those of 
the nutations but whose coefficients change linearly with time. Regard ing the 
nutation, the following observations are made from (2.2-42) and (2.2-43): 


1. Each nutation constituent is in general an elliptic motion. 
It can be separated into two circular motions which have 
opposite angular velocities . If Aj- - Aj + , the nutation is 
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circular. There- will Be no term in the obliquity. An 
example of' such a 1 motion is the annual nutation term 
due to the solar tides. 


2. Because of the presence of A aj in the denominator the 
contribution of a particular tidal wave to the nutation 
depends not only on the magnitude of the tidal wave but 
also on the location of its frequency relative to the sidereal 
frequency. Those waves which are close to the sidereal 
frequency 165.555 have an amplified effect on the nuta- 
tions. 

3. The nutations are officially tabulated in the form 


x HN < aj cos Aa -i+ + ib J sin Aa J+ ) 

j 

A j + -Aj _ 


(2. 2-46) 


aj = - 


with 


bj = 


A3 £2 A a j + 

1 Aj + + Aj _ 

A3 £2 A dj + 


It is customary in astronomy to count the Euler angles 
ip positive westward. This is contrary to the convention 
adopted here. 


4. The method chosen here for deriving the nutation, i. e. , 

> » 

using the precise numerical expressions from Doodson's 
development, concedes an important relationship between 
the constant of nutation N and the luni-solar precession 
ip H p. N is the nutation coefficient in obliquity associated 
with a 18.6 year period. Using a series expansion for 

' ^ " .. * w - ' 1 

the eccentricity and inclination when, computing the 
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lunar and solar coordinates the following literal expres- 
sions can be shown to hold [Kulikov, 1956]: 


A3 -Ai 

^hp = 7 cos 6 Hp (K f~ t + K') 


A3 
A3 - Ai 


l+/i 


— (4871'.' 140 + 866195V 0 • -^-) 

A3 b-fx 


N 


H cos 6 


HP 


1 +/2 


A 3 -Aj 
A3 


A3 — Aj 


231 981.8 


._M_ 

l+/i 


(2.2-47) 


I u is the mass of the moon in units of the earth mass and 
the coefficients K, K 7 , and H ' are functions of the orbit 
elements of the moon and the sun. The latter can be com- 
puted from theory with sufficient accuracy. The numer- 
ical values which are given here refer to the standard 
epoch 1900.0. Combining equations (2.2-47) gives 


H 


N = 


tf>HP ( 




K 


K 

47.6237 f i 

1+ 178.822 
r 


fl K 


■f-1) 


-1 


(2.2-48) 


The importance of these formulas is that given the observa- 
tions for i/J H p and N one can compute the dynamical flatten- 
ing and the mass ratio ix. However, the mass ratio obtain- 
ed in this manner is not consistent with the values based 
on recent data from lunar and planetary space craft. This 
discrepancy is one of the outstanding problems in the sys- 
tem of astronomical constants. It indicates that the under- 
lying ample nutation theory is not adequate. A partial answer 
is given when considering earth models with liquid core 
(See Section 3). 
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Relative Motions 


2.3.1 The Complete Excitation Function 

In order to proceed with the solution of equation (2.1-9) for the 
motion of the rotation axis (I) the disturbances c in the inertia tensor 
and the relative angular momentum h, which are to be used in the excita- 
tion function (^) need to be computed. For a completely elastic body 
(Hooke body) the derivations can be found in the standard literature. The 
following expressions for the disturbances of the inertia tensor are taken 
from McClure (1973). 

The perturbations due to rotational deformation are 
k 

C 13 rd = ^ (A3 - Ai) U U 


c 33 R 0 


k 

k e 


(A 3 - A l ) u 3 i 


In complex notation the perturbation is written as 
k 

Cro = ^ (As - Aj) ui 

The secular Love number k s is 

= 3G(A 3 - A : ) 

a 5 a 2 


(2.3-1) 


(2.3-2) 


and k is the tidal effective Love number. The perturbations due to tidal 
deformation are 

kM D 

C 13 to = a 5 ^ A 3j j sin a } 


C 33 Tp 



a 5 ]TA 31J 

j. 


cos a j 


or 
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(2.3-3) 


C T 0 


kM, 


= — D - a^lAaue 


-Wj 


j 


The misalignment of the third axis of (U) with the maximum moment of 
inertia axis (F) is expressed by the constant corrections 


Cq — c 130 "** 1 C330 


(2.3-4) 


Only those perturbations are given above which directly effect the 
equations of motions (2.1-9). These perturbations result from the second 
degree tesseral harmonic coefficients of the expansion of the earth' s 
gravitational and tidal potential. The zonal coefficients effect only the 
diagonal terms of the inertia tensor, and according to equation (2.1-10), 
the velocity of rotation of the earth. The sectorial coefficients contri- 
bute only c n , Cg, , and c 12 » The excitation due to the relative angular 
momentum h in is according to equation (2. 1-9) 



h 

(A 3 - Aj)fl 


ih 

(A 3 -Ai)f2 2 


An estimation of ^ m the case of an elastic body is given in [Heitz, 
1976, Appendix 3] . He found that ^ is negligibly small. 


2.3.2 B ody- Fixe d Motion 

The complete solution for the body-fixed motion of the instantane- 
ous rotation axis is the sum of the homogeneous solution (force free, 

'& = 0) and a particular solution (forced solution, 'l' ^ 0) of the differen- 
tial equation (2.1-9) 

U! = i CT r (ui - St ) 

The complete solution is 

t 

icr r t icr r t r -r -ia r T 
uj = uo e - e iCT r / ^ (T) e dT (2.3-5) 
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u 0 is a complex constant of tHe homogeneous solution. The excitation 
function & in (2.1-9) is given by equations (2.2-20), (2.2-28), (2.2-29), and 
(2.3-1) to (2.3-4). The solution for Ui after integrating (2.3-5) and 
some lengthy algebraic rearrangement is 


ui 


u 0 e 


i(7 jjt 


* b + 1 


E 

J 


s' 1 . A.e 
A- 


-l a , 


(2.3-6) 


where 


‘ (To 






i 




a 

a 


i + 


Aa 

n 


k 

k. 


+ 


av 

0 


(2.3-7) 


The motion of the angular momentum axis (H) in the system (U) is 
according to equation (2.1-12) 


u n 


Ar 

A 3 


Ut + 



+ 


h 

a 3 a 


(2.3-8) 


The term due to the relative angular momentum h can be neglected again. 
Thus, substituting the expressions (2.3-6), (2.3-1), (2.3.3) and (2.3-4) m 
(2.3-8) gives the body -fixed motion of (H) 


where 


U H 








1 - 
a 


k 




% + 


E 


s '< s j 




(2.3-9) 


(2.3-10) 
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The motion of the axis of figure is according to equation (2.1-16) and equa- 
tions (2.3-1) - (2.3-4) 


u 


F • 


k 

k. 


u D e i<Tot + 


ik y sf B " A 

J P 3 (A a - A 1 ) J 


(2. '3-11) 


where 



(2.3-12) 


The basic equations (2.3-6), (2.3-9), and (2.3-11) for the respective 
motions of the instantaneous rotation axis (I) , the angular momentum 
vector (H), and the axis of figure (F) can be re-written in terms of the 
nutation argument Aoq with the help of equation (2.2-32). The forced 
components are 


u 


if 



s‘f -i (GMST + Aa 3 ) 

Aj 0 

a 1 o 3 


u Hf E 


s T s, 


'S i A o 


-i(GMST + Aa^) 


_ * i5 

u "~\ E 


-1 " 
S J S 3 


A. e 


-i (GMST + A o^) 


3 


, 0 3 (A 3 -Ax) i 


(2.3-13) 


The homogeneous, or force-free, motions have a conceptual import- 
ance in defining polar motion. They give the hypothetical positions of the 
axes which would be occupied if there were no external forces acting on 
the earth. These motions, already contained in the complete solutions, 
are separately given by 
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Uih 

UHh 

l 

Ufh 


7o e 

r 

yo 


i(a o t + r ) + ^ 


Ai 

As 


A 3 - A i\ k 
A 3 


i(ff 0 t + F) + ^ 


(2.3-14) 


k i(a o t + r ) , T 
7o “ e v ’ + $0 

k s 


where the complex constant u 0 is written as 

iT 

u 0 = Yo e 

F denotes the phase angle for t = 0. 

The direction cosines u are transformed to the local components, 
defined by the direction of the local meridian and the direction orthogonal 
to it, by 

♦ 1 

u' = u e 1 (2.3-15) 

An equivalent procedure is to replace GMST by mean sidereal time MST 
in equations (2.3-13). Then, the along-meridian component is Re(u ) ~ 
and the across-meridian component is lm (u 7 ) = . 

Figure 2. 3 displays the body-fixed motions for the elastic model 
as discussed so far. The figure represents a tangent plane on the unit 
sphere, the point of tangency is the point where the third axis (U 3 ) of 
the coordinate system (U) passes through the sphere. Each of the 
motions indicated is the equivalent of a mathematical expression 
developed above'. The following observations can be made: 

a) The separation U3 - % j which is given by the second equation of 
(2.3-7), is a linear function of the constant component c 0 in the 
inertia tensor. (#□) is not identical with the axis of figure (F). 

For the rigid (k = 0) model the axis of figure (F) will be at T 
at any time. 
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b) The radius of the circle with center at 4> 0 and which passes through 
(S) is given by the first term of u Ftl in equation (2.3-14). Thus, 
the position of (S) is given by u f h in (2.3-14). Similarly the 
points (C / ) and (E 0 ) are found by plotting u dh and ui h , respectively. 

E 0 has been termed in Woolard [1952, p. 160] the "Eulerian Pole 
of Rotation" . It corresponds to the hypothetical surface point which 
the instantaneous rotation axis (I) would occupy if no external forces 
were present. The prograde (in the sense of the earth rotation) 
body -fixed motion of (E 0 ) and (C / ) is sometimes called "wobble". 

The amplitude and phase cannot be predicted from theory. The 
period is either the Chandler (elastic) or the Euler (rigid) period, 
as is well known. 

c) The nearly diurnal motion of the axes (I) and (H) around (E 0 ) and (C) 
are retrograde motions. They are a result of luni-solar attraction 
as expressed by equations (2.3-13). In case of the elastic body even 
the axis of figure (F) has a non-zero diurnal motion u Ff . The 
motions m Figure 2.3 are approximated by the exact circular diurnal 
constituent corresponding to the tide Kj which causes the precession. 
In actuality each tidal frequency has a corresponding diurnal motion 
component. Numerical values for these motions are given m McClure 
[1973] and Woolard [1952]. The equations (2.2-29) and (2.3-13) show 
that these forced motions are proportional to the difference in the 
principal moments of inertia. 

d) The complete solution as well as the homogeneous and forced solu- 
tion separately fulfill the relation (2.1-19) 

Ai , > 

U H ~ Up = -J— ( UI - U F ) 

as can be verified from the respective equations. Thus, all three 
axes are in a plane at all times. 

The axes through the points (S), (C'), and (E 0 ), which result from 
the homogeneous solution, exhibit no periodic diurnal body-fixed 
motions . 40 
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Figure 2. 3 Body- Fixed Motions for the Elastic Body 
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e) The magnitudes of the various motions in Figure 2.3 are given in 
Table 2. 1. 


Table 2.1 

Magnitudes in Polar Motion 


Forced Motion 

Model 

I - Eo 

I - H 

S - F 

elastic 

^ 62 cm 

^ 21 cm 

^ 60 m 


rigid 

^ 61 cm 

1.5 cm 

F = T 


Free Motion 

Model 

Radius 

Frequency 

C 1 - Eo 

elastic 

7o 

(To 

( l 

W 

rigid 

7o 

ov 




f) The separation E 0 - C' is due to the difference in the homogeneous 
solutions of the rotation axis (I) and the angular momentum axis 
(H) . Its magnitude 


7o 


'Ab-aA 
, A, / 



=s 1.5 cm 


can only be indirectly obtained through the observed value of y 0 • 
In the next section this distance will be verified as the amplitude 
of the free nutation. 


As already mentioned the -motions of the rigid earth model are obtained by 
putting k = 0. Equations (2.3-6), (2.3-9), and (2.3-11) are then 
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Ui. 


Hk 


Uf 


i(o v t+r) Y' ^ e" ia - 

Ai i(cr r t + 37) 
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yo~ e 
A 3 


* = , .- C ° , , 
r A S - Ai 


+ 


-*£ 




-i a , 


As 0(a j + <x r ) 


(2.3-16) 


For a rigid model the axis of figure (F) is, of course, body-fixed. Figure 2.3 
can easily be modified for the rigid body ease. The separation between the 
instantaneous rotation axis (I) and the angular momentum axis (H) is much 
smaller than for the elastic body because the coefficient s ' s in u rif of equation 
(2.3-13) is equal to 1. The following' ratio holds: 

u^. , elastic _ _ k 

u Hf , rigid k s 

The "forced" terms of uj ( are usually called the Oppolzer terms. The 
most important ones are for the rigid model [Woolard, 1952]: 


Unf “ 0.0087 sm (GMST) 

-0.0062 sin (GMST-2s) 
-0.0029 Sin (GMST-2h) 
+0.0012 sin (GMST-O) 

vi 2 if “ 0.0087 cos (GMST) 

-0.0062 cos (GMST-2s) 
-0.0029 cos (GMST-2h) 
+0.0012 cos (GMST-fl) 

The arguments are explamed in (2.2-22). 


(2.3-17) 


2.3.3 Space- Fixed Motion 

The spatial positions of the various axes are determined with respect 
to the angular momentum axis (H) , whose position ( 0 H } ) is known 

from the solution of Poisson's equations. The direction cosines uj. which 
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were derived in the previous section are transformed into changes of the 
corresponding Euler angles. 

Consider the previously used coordinate system (X) H and (U). The 
position of the first system relative to the system (X) E is given by the 
pair of Euler angles (0 H > 0 n ) and the relation (2.2-9) 

(X), = Rs H h ) Ri (S„) (X)* 

The body-ftxed system (U) is related to the inertial system by 

(f) e = Ra H ) R a ( 9 ) Ra (-<p) (U) 

These two equations readily give the direction cosines of the axis (H) in 
the (U) system. For small differences in the Euler angles the direction 
cosines of the angular momentum vector (it) are to the first order of 
small quantities 

u 1H 6 ^ sin 0 h 

u^h = Ik{<p) -60 (2.3-18) 

U3 H J L 1 

where 

bij) = ij) - 
60 = 0 - 0h 

The first two equations of (2.3-18) can be written in complex notation as 
follows : 

u a H + iUa H - e l(fi (6 ij) sin 0 K ~ 16 0) , 
which is rewritten as 

= 60 + i 5ii>sin0-, - i e ic0 (u lH + iUg H ) (2.3-19) 
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Note that the correction is given in the .sense of 




XH 


where the subscript U 3 is used to emphasize that x^ determines the 
space-fixed position of the U 3 - axis relative to H. 

The space-fixed position of the various axes can now be given. 

Using (2.3-19) the instantaneous rotation axis (I) is obtained as follows: 

Xj = x H + 6 x T 

= X H + 6x Ih + Sx If 

- X H + 6% + i 6 il)r sin 0H 
= x H + ie *P(u H - Uj) 


Substituting equations (2.3-6) and (2.3-9) the Euler angles become 

-i Aaj 




(2.3-20) 


where 


B, 


fA 3 A, 1 _ j 


A n 


s j (A 3 -A 1 s') 3 

A^O 8 


The space-fixed motion of the axis of figure (F) is with equations (2.3-9) 
and (2.3-11) 


=x H + 5 x f 


. Ai / k\ i( CTo t+r-Hp) v- -iA“j 

+ ? B ’ e 


(2.3-21) 


Finally, the spatial motion of U 3 - axis is 

x =x + 5x 
U3 H U3 

= x H + ie^ u H 


= x H+ i 


L yo e l(g ° t+ 1 >(p) + y^ S A 3 e lA0!j (2.3-22) 
II j A3O 2 


A^ + \ A 3 


+ i e' 


i<P 
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Note that during the derivation of these equations we have used the 
approximation (2.2-37), i.e., the Euler angle <p has .been equated to 
Greenwich sidereal time. 

The following observations can be made regarding the space-fixed 
motions: 

4 

a) The homogeneous term 

6x rr = -j ^ A 3 ~ _ k j y ne i(aotTli+(P) 

of equation (2.3-20) is called free nutation . It gives the hypo- 
thetical path of the instantaneous rotation axis (I) in space which 
would occur if no external forces were present. Its frequency 
is nearly diurnal. The amplitude, based on the observed values of 
y 0 , is approximately 1.5 cm. The free nutation is not included 
in the officially adopted set of nutations of the rotation axis . 

Special considerations will be given in Section 2.4. 2 as to the 
observational significance of the free nutation. In any case, a frame 
connected to the instantaneous rotation axis (I) will perform a nearly 
diurnal "rocking 11 in inertial space. 


b) The forced periodic terms in equation (2.3-20) which are to be added 
to those of equation (2.2-43) of the angular momentum axis (H) in 
order to obtain the forced motion of the instantaneous axis of rotation 
(I), are of significantly large magnitude for the elastic model. They 
account for the 21 cm separation between (I) and (H) already indicated 
in Figure 2.3. The ratio with respect to the motion of the angular 
momentum axis 


16 x 


r? 


1 


x. 


A3 1 

v 1 n } 

l£ 

~ Ai| 

(l k j 
v 1 k J_ 


Ai 

aj + 

C r ( 

h + k \ 
\ TTiZJ 

- 


demonstrates the dependencies on the elastic properties. 
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e) The axis of figure (F) has, of course, a diurnal nutation component 
too, which results from the homogeneous solution. The forced 
motion is very large; at describes the 60 m motion indicated in 
Figure 2.3. The axis of figure is the most sensitive axis with 
respect to mass redistribution among all the axes considered so 
far. It is, therefore, unlikely to be chosen as a defining direction 
for a reference frame. 

d) The equations give immediately the ratios: 


6x Ih 

= 6*1, = A, - A, 

6 Xf}i 

Sx Ff 


The motions for a rigid model are obtained by putting k = 0. They are: 




In this case the ratio 

. _ A3 - Ai ^ Aq?j 

x h ^ A3 a, + a r 


is extremely small. The largest periodic term in Sx# has a coefficient 
of 0”00002 .[Woolard,, 1952, p. 133]. The terms are, consequently, not 
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included in the officially adopted set of nutations for the instantaneous 
rotation axis (I) since terms smaller than. 0'-' 0002 are omitted [AENA 
Supplement, 1961, p. 44]. The coefficients Bj for the rigid model are 
given in Woolard [1952, p. 132]. The largest terms are 


0.01615 

sin(2s) 

+0.00753 

sin( 2h) 

+ 0.00338 

+ ... 

sin (£1) 

-0.00868 

+0.00590 

cos (2s) 

+0.00275 

cos (2h) 

-0.00100 

cos (O) 


(2.3-24) 


The arguments are explained in equations (2.2-22). The constant term in 
5 0^ f is caused by the sidereal term in 6* ff , it is actually a combina- 
tion of two sidereal terms, one resulting from the moon, the other from 
the sun. The motions 6x rf and are referred to as "luni-solar 

diurnal nutations" in the astronomical literature [Woolard, 1952]. This 
terminology seems an unfortunate choice since these motions have no 
diurnal period whatsoever. 


2.3.4 The Celestial Pole (C' ) 


The pole to which the nutations refer is denoted by (C ') . This 
pole fulfills the following criteria: 

1) The pole (C / ) should exhibit no nearly diurnal periodic body-fixed 
motions . 

2) The position of the (C') in space should be computable at any time 
from the motion theory of the underlying models, i.e., no free 
solution component is permissible. 

"For ideal bodies, such as the rigid or perfectly elastic earth 
model-, both criteria are easy to fulfill. Considering Figure 2.3, there 
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The expressions are arrived at by identifying the body-fixed components 
in Figure 2. 3 and using their transformed values in terms of corrections 
to the Euler angles. 

Only the celestial pole (C / ) qualifies under criterion 2. Its position 
in space can be computed from equation (2.2-43) for X H and from equa- 
tion (2.3-22) for 6x Usf . The north celestial pole thus defined has no 
periodic diur nal body -fixed motion due to external forces. Its body-fixed 
motion results solely from force-free motions. Its space-fixed motion 
is entirely due to external forces; it has no space-fixed motion resulting 
from a homogeneous solution component. The point, E 0 , which is called 
by Woolard the Euler ian pole of rotation, exhibits a diurnal motion in 
space which is exactly equal to the free nutation given in equation 
(2.3-20). The amplitude and phase of 6x rb cannot be determined from 
theory. 

The term 6x Jaf is the transformation of uu f in equation (2.3^13). 
As equations (2.2-29) and (2.3-22) show, the amplitude of 6x Ugf is pro- 
portional to the mass of the disturbing body Mjj and the form factor 
(A3 — Aj )/ A3 . Another interesting feature of the pole (C') is seen if we. 
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for reasons of simplicity, assume a zero free motion. This is the case 
of the steady state motion. The pole (C*) has no body-fixed motion in 
this case. But, the angular momentum axis (H) and the instantaneous 
rotation axis (I) still have the nearly diurnal motions u hf and u rf , 
whose magnitude is proportional to and (A 3 -Aj )/As . If one constructs 
various models with different form factors, i.e. , adding mass points 
symmetrically around the equator, the pole ((/ ) will remain at the same 
body -fixed position whereas (H) and (I) change the radius of their circular 
motions. It is thus clear that the celestial pole {d ) and the instantaneous 
rotation axis (I) do have different properties and should conceptually always 
be separated. However, from a merely descriptive point of view, the 
pole (C') may be considered as a ’’rotation axis" as well. 

It is worth noting that in the case of the rigid earth model the 
spatial motion of the north celestial pole is 


K,r . . , = + 6tt F , . . , 

c rigid H F rigid 

since according to equations (2. 3-23) the equality 

6 x - 6 X Ff ..j 

u 3 frigid Ff rigid 


(2.3-26) 


holds for the rigid body. The spatial separation between the pole (C 7 ) 
of the elastic and the rigid model is 

A(V) = 6 *t' 3 f elastic " 6 \f rigid 

The maximum separation is 21 cm (see Figure 2,3). 

It is probably because of equation (2.3-26) that in the literature one 
finds frequently the statement that the nutation should be computed for 
the axis of figure [Atkinson, 1973, 1975, 1976]. 

There is no doubt that in all 'eases the axis as expressed in equa- 
tion (2.3-26) was intended. But, once again, for clarity’s sake, one 
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ought to distinguish between the axis of figure (F), which is defined as 
the direction of the maximum moment of inertia at any epoch, and the 
north celestial pole (C / ) as has been defined above. Even for the rigid 
earth model the axis (C ) generally does not coincide with the axis of 
figure (F). Only when the homogeneous solution is zero the identity 

x F = * c , 

holds. In case of rigid motion the intermediary role of the instantaneous 
axis of rotation becomes quite obvious. When setting up the Euler dynamic 
equations (2.1-3) or (2. 2-5)the concept of instantaneous rotation axis is 
needed since (A are the velocity components of the instantaneous rotation 
axis (I) in the frame (U)f. Symbolically Euler’s dynamical equations can 
be written as 

EDEq= fi(w 1 .lr) = L if 

Euler’s geometric equations relate the velocity eomponets co t to the 
derivatives of the Euler angles 0 F ■ , >b~ ,, <jy. The Euler angles relate 
the frame (U)f to the inertial system. Substituting 

V 

• • 

Wj — Si ( 0f>^>f 5£ 0f) 

in the dynamical equations gives 

EDEq = f t (g (Bp , jJ) F , <£,>,!,),= L lP (2.3-27) 

These equations can be solved for the Euler angles.' The angles 0 f and 
ip F give the ’’forced position" of the axis of figure in space. For a rigid 
body the axis of figure (F) has no body-fixed motion, hence 

Up ? rigid 0 

The space-fixed motion is strictly computable from equation (2.3-27). 
Therefore, the so-defined axis fulfills both requirements for the pole (C') 
which were set up above. In fact, the solution of equation (2.3-27) is 
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identical to equation (2.3-26). One can now add the homogeneous solution, 
which results in the deviation of the ((/) from 'the axis of" figure, but this 
aspect is not of concern here. Eather, it is noted that after substituting 
Euler’s kinematic equation into equation (2.3-27) one does not have to 
worry about the instantaneous rotation axis (I) anymore. The instantaneous 
rotation axis (I) is conceptually needed only at the initial stage when 
formulating the dynamical equations. 

At present the rotation axis (I) is adopted as the reference pole 
instead of (C'). The practical complications arising from that convention 
are dealt with in Section 2.4.4: In Section. 2.4.% the observational signifi- 

cance of the homogeneous component which separates ((/) and (Eo) in 
Figure 2.-3 is discussed. The additional complications arising from the 
liquid core are the subject of Section 3. 

2.3.5 Poinsot’s Kinematical Representation 

Poinsot showed in 1857 that a continuous rotational motion of a 
rigid body about a fixed point is always geometrically equivalent to the 
rolling without slippage of a body cone (polhode cone) on_ a space-fixed 
cone (herpolhode cone). The line of contact between these two cones is 
the instantaneous rotation axis (I) of the body. Mathematically, Poinsot’s 
representation is related to the motion of the inertia ellipsoid. An 
extensive treatment on this subject is given in the standard literature. 

Here only a simple intuitive explanation is given. 

Any .continuous motion of a body about a fixed point can be repre- 
sented by dividing the time into infinitesimal elements and considering 
the motion of the body during each element of time as a rotation about 
the corresponding instantaneous axis of rotation. During the motion, the 
position of the instantaneous axis in the body and in space is gradually 
changing. To show this, consider the case of a cone rolling on a plane 
and rotating about its fixed vertex 0 (Figure 2.4); 
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Figure 2.4 Poinsot’s Kinematical Representation 


At any instant, the motion of the cone consists of a rotation about 
the line of tangency OA, which is the instantaneous axis in this case. As 
the cone rolls, various lines of its surface come into contact with the 
plane , and the instantaneous axis has different positions in the moving 
cone. At the same time, consecutive positions of the instantaneous axis 
in space form the plane on which the cone is rolling. This example can 
be generalized. Instead of a circular cone rolling on a plane, a cone of 
arbitrary shape is rolled on a surface of another cone. One can assume 
also that the rolling cone is not a physical body but a geometrical sur- 
face formed by consecutive positions of the instantaneous axis in the mov- 
ing body of any shape as indicated by dotted lines. Likewise, the space 
cone is formed by successive positions of the instantaneous axis in space, 
and the motion of the body is visualized by rolling the cone connected 
with the body on the cone fixed in space. By varying the shapes of both 
cones, all possible motions of a rigid body about a fixed point can be 
obtained. 
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In case of the forced motion of the earth, the Z-axis coincides with 
the pole of the ecliptic and the space fixed cone has a mean radius equal 

to the obliquity .. The cone is not circular but has -"ripples" superimposed 

< 

on it due to the nutation. The time varying size of the body-fixed cone 
is determined at any instant through the diurnal polar motion terms of the 
rotation axis, e.g., ut? of equation (2.3-13): It is understood that the 

term "body-fixed" cone refers only to an infinitesimal time span. With 
this in mind, the two cones describe the forced position of the instantane- 
ous rotation axis completely at any instant. The center of the body cone 
can be associated with the pole (C'). But it is known that the body cone 
or its representative, the axis (C'), moves within the earth due to the 
Chandler motion (force-free motion). This motion itself can be represent- 
ed by a pair of rolling cones. This leads to an alternative representa- 
tion of the motions, i.e., instead of using one pair of general cones in 
infinitesimal time intervals several pairs of circular cones are used in a 
finite time interval. This allows the separate interpretation of the various 
motion components. Each of these pairs of cones defines a particular 
component cOj of the instantaneous rotation axis in magnitude and in 
direction. The real physical instantaneous rotation axis in space is the 
sum of all components: 


~ 5 > 

CO 


E 

j 




Besides separating the free and forced motions the latter is once more 
split up in precession and circular nutations. Thus 

w = N 

3 

It is more convenient to deal in this context with the circular nutations 
rather than with their elliptical combination since the former are directly 
accessible to the cone representation. As a type of classification a 
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motion is called prograde if it occurs in the direction of the earth rotation, 
otherwise it is called retrograde. The three types of motions which occur 
in the earth motion spectrum are shown in Figure 2.5. The angles s and 
r are the vertex semi-angles of the space-fixed and body-fixed cones, 
and v and p are the angular velocities of the instantaneous rotation axis 
component on the respective cone. The vertex of all cones is at the 
geocenter. The magnitude of the resulting rotation axis component is 

u> i = p 2 J v* - 2 [iv cos 6 

i 

The two equations, as depicted in Figure 2.5, 

* * ■» 

sin r _ sin 6 
v u; 

sin s sin 6 
p (jo 

give the relation between the vertex semi-angles and the angular velocities 


sin r _v_ 

sin s p. 


(2.3-28) 


This equation implies that the cones move on each other without slipping. 
For small vertex angles the relation simplifies to 

r v 

- = - . (2.3-29) 

s p 


The various motions can now be classified. 
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A) Force-free Motion: Case 3 


Z = angular momentum axis (H) 

N = axis of figure (F) 

From the first equation of (2.3-14) and equation (2.3-20) the 
ratio (2.3-29) is Vq a 0 + i Q 


V 0 




CTo 


(2.3-30) 


Here we have used Q,= <p which results from the approxima- 
tion GMST = (p. The validity of (2.3-30) can be verified by 
rearrangements on the left hand side. 

B) Precession: Case 1 ; = 0 

Z = fixed ecliptic pole 
N - c' (precession only) 

From uif of equation (2. 3-13) the vertex semi-angle of the body 
cone for A a* = 0 is 

i 



— 


The vertex semi-angle of the space cone is 0 O as is seen from sim- 
ple geometric considerations. The spatial frequency y is taken from 
equation (2.2-44) as 

; _ A k 

V ~ 0„P " 11 

A 3 £2 sin 0 Q 

The, body-fixed frequency is, of course, equal to the sidereal fre- 
quency O. Thus, equation ,(2.3-28) gives 

Ak x ^ 

^ 2 A 3 A s ^sin0 o 

sin 0o 
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C) 


Nutation: Case 1 if AcE^ >0 

Case 2 if Aa^ < 0 

Z — is located on the space cone of the precession 
N = C f (nutation only) 

Attj is the nutation frequency. It was first defined in equation 
(2.2-32). The nutation frequency is either positive or negative, de- 
pending whether the corresponding tidal frequency is larger or smaller 
than the sidereal frequency n. The vertex semi-angle and frequency 
for the hody-fixed motion is given by equation (2.3-6), and those for 


the space-fixed motion are taken from, equations (2.2-43) and (2.3-20). 
Inserting these quantities in (2.3-29) gives the ratio: 


A i s ? 

Al _ Aoc i (2.3-31) 

~T S A J s“ x (A 3 - A x s') flij 

, -f 

A 3 QAa i Q 3 A 1 A 3 

The validity of equation (2.3-31) can be verified by re-arrangements 
on the left-hand side. The equation holds for the rigid and elastic 
models. 


In general, then, all circular nutations and precession will result in 
retrograde body motions because the exponent in the expressions for u If 
im equation (2.3-6) is always negative. Similarly, those circular nutations for 
which Ac^ < o will cause prograde spatial motions of the rotation axis, 
whereas the others give retrograde spatial motions. 

It has been demonstrated above that all known motions for the rigid 
and elastic model can be ‘represented by pairs of cones rolling on each 
other without slipping. Each type 'of motion is completely identifiable by 
its period and its occurrence in either space or body. 
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2.4 Coordinate System Definition and Observability 


2.4.1 Relationship Between Coordinate Systems 
The Poisson equation for the spatial motion of the angular momentum 
axis was solved in an inertial frame which was taken to be the ecliptic at a 
fixed (standard) epoch. The physical realization of the ecliptic can be thought 
of as the plane defined by the solar center, the earth-moon barycenter and the 
velocity of the barycenter. The orientation of this plane changes secularly and 
periodically due to the planetary perturbations. The periodic variations are 
computable from theory. That fictitious plane which has only a secular motion 
in inertial space is the exact definition of the ecliptic. 


The relationship between the ecliptic and the celestial system whose 
pole is (C') is given by the dynamical theory of the rotation of the earth 
in combination with some constants of definition. 

The mean position of the celestial equator (C^) relative to the fixed 
ecliptic is implicit in equation (2.3-25). Using equation (2.2-45) and the 
constant term of dx^* in (2.3-22), which is identical to the constant 
term in (2.3-24) in case of the rigid body, the precession of the celes- 
tial pole (C / ) is 


0c'p = 0 hp = fit + f 2 i? + . . . 

0 C ' P = 9hp + 60° = 0o + 66 * + 0 2 1 2 + . . • 


(2.4-1) 


f x is the constant of luni-solar precession. It can be derived from ob- 
servations or from theory. The latter possibility is indicated through equation 
(2.2-47). The mean obliquity at the standard epoch is found by obser- 
vations. This constant includes the small constant d0 c of equation (2.3-22) 
or (2.3-24) which is equal to the mean separation between the angular 
momentum (H) and the celestial pole (C 7 ). For abbreviation, the equations 
(2.4-1) can be rewritten as follows: 

X(pp — Xhp + 6 X°u • (2.4-2) 

3 

The superscript c denotes the constant term in 6 Xu t • 

3 
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The nodal line T , in Figure 2.6 between the fixed ecliptic and the 
mean equator of date does not coincide with the direction of the mean-vernal 
equinox T M , The latter is equal to the nodal line of the ecliptic of date 
and the mean equator of date. The motion of the ecliptic on the moveable 
mean equator is called the planetary precession in right ascension a, with 

a = g x t + g 3 t 8 + ... 

The g x 's follow from the planetary theory. The transformation due to 
precession is 

s ~ ^3( a ) ^(-©c 7 ?) H>c' P ) Hxtflb + 5 0°) {2.4-3) 



and the remaining transformation from the mean celestial system to the 
true celestial system <X) C/ is 


where 


N = B 3 (a 1 ) R3(-0 C / P - 

(2.4-4) 


a x = a + A a 

and Aa is a small computable term [Woolard and Clemence, 1966, p.240]. 
Gc'n and 3X6 the nutation angles of the celestial pole (C y ). With equa- 
tion (2.3-25), the nutation angles can be written as 

*c'n " *hn + 6 Hoi t 

3 
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where the superscript denotes the periodic terms in 6 Xu t . 

3 

With a few minor changes. Figure 2,6 can also serve for setting up equa- 
tion (2.4-4). 

hi practical astronomical work one tends to select a different set 
of rotation angles in the transformation (2.4-3) [Mueller, 1969', p. 65). 

One also prefers to use the so-called 1 'reduced" nutation in longitude 
and obliquity, ^)'/ N and 9 q/ n , which relate the true equator to the 
mean equator of date. Then the transformation N becomes 
N = ~ 

The reduction of the nutations is carried out in Woolard (1952, p. 157). 
e is the mean obliquity of date, 

e = 9 0 + 8 0 C + a-jit + a 2 t 3 + . . . 

The coefficients a i follow again from theory . 

Finally, the motion of the vernal equinox due to the combined motion 
of the ecliptic and the equator is called general precession in longitude 
ibi> 

0! = h x t + h^ t 2 + ... 

where 

h x = f x - g x cos(0 o + 6 0 C ) 

Since g t and h 1 follow from theory there is no additional new constant 
which is to be determined by observation. 

2.4.2 Fundamental Declinations and Latitudes 

The fundamental declinations and latitudes are the true observables 
in observatory astronomy. They give the position of the observatory 
and the star in the true celestial system (X) c /. The procedure consists 
of observing zenith distances of stars at upper and lower culmination. 

At no stage of the procedure will the star position or an adopted series 
of nutations be needed. 
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For the following description it is assumed that the Chandler motion 
is zero for the period of 12 hours. This assumption, of course, will 
not be fulfilled since one has no control o.v.er.the actual free motion. 

This simply demonstrates that all efforts to determine the orienta- 
tion of the earth in space are ultimately limited by the amount of Chand- 
ler motion which occurs during the interval which is needed to carry 
out the basic observations. The progressive Chandler motion is approx- 
imately 5 cm for 12 hours. 

The left and the right picture in Figure 2.7 show the situation 
at upper and lower culmination of the same star. The pole (C') remains 
body-fixed during the 12 hour interval. In particular, it has no nearly 
diurnal free nutation as was shown in Section 2.3.4. The body-fixed 
motion of the angular momentum axis (H) is represented as a circle in 
the lower two pictures. There is no need for (H) or the axis (U 3 ) to 
be in the meridian during the observations. The body -fixed position 
of (H) for upper transit at epoch T is shown in the lower left picture. 

The lower right picture shows the body-fixed position of (H) at lower 
and upper transit. For the sake of completion, a possible body-fixed 
position of the Eulerian pole of rotation (E 0 ) is also shown. The 
angle y is the along- meridian polar motion component . 7i and y 3 are 
the periodic diurnal body-fixed motion components of (H) along the 
meridian. ^ denotes the adopted station latitude in the (U)-system. 

6 is the declination in the true celestial system (X) c / whose third 
axis is identical with ( C '). The observed zenith distances at culmina- 
tion are z z and z Q . From the upper two pictures in Figure 2.7, the 
following two basic relations are readily seen: 

90-6 = <z 2 - Zj)/2 


(2.4-5) 


90 - ($u 3 + y) = (z a + z^/2 
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It is now clear that the declination and the latitude in the true celestial system 
(X)c'are the observables . The position of (H) at the time of observation has 
no importance in this context. More importantly, the instantaneous rotation 
axis (I) does not enter the observation geometry at all. It is therefore concluded 
that, on the basis of the observational procedure* discussed above, the instan- 
taneous rotation axis (I) and the angular momentum axis-(H) are not observable. 

The observational insignificance of the Eulerian pole of rotation 
E 0 is readily demonstrated with the help of Figure 2.7. As has been 
pointed out previously, the only difference between the north celestial 
pole (C*) and the Eulerian pole of rotation (E 0 ) is a homogeneous solution 
component, so that the latter still exhibits a nearly diurnal space-fixed 
motion due to free nutation. Consider the body-fixed point E 0 in Figure 
2.7, In the upper figures, E 0 performs a similar motion relative 
to the ( C *) as does, for example, the axis (U 3 ). But the presence 
of such a fixed crust point, and there are infinitely many possible 
choices all laying on the line \& 0 - C' of Figure 2.3, does not effect 
the principles expressed by equations (2.4-5) at all. Thus, the Eulerian 
pole of rotation Is not observable either. 

2.4.3 Determination of the Constants of Definition 

The relationship between the celestial system (X) c / and the fixed 

ecliptic system (X) e is expressed by the standard expressions 

sin 6 = cos ft sin> sin0 c /f + sin# cos Q c f 
and ( (2.4-6) 

cos 6 sin O' = cos R sin \ cos 9 C / - sin g sin0 c / 

and 6 are the right ascension and declination in (X) c , . \ and R are 

the ecliptic longitude and latitude in the fixed ecliptic system. Q./ is 
the true obliquity of date, which is split up according to equations 
t2.4-l) and (2.4-2) as 

Qc 1 ~ 6c' p + 9c' n (2.4-7) 
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whereby the small constant 6 0 C is included in Q c , p . Eliminating the 
longitude \ in (2.4-6) gives 

sinatangj., = tan 6 - Sin ^ (2.4-8) 

• . cos 6 cos Q c , 

The ecliptic latitude g of objects close to the ecliptic such as the sun, 
planets, etc. are known from theory accurately enough m order to be 
treated as a known quantity. Since at the meridian transit the right 
ascension c t is equal to the sidereal time, the right ascension is ex- 
pressed in terms of the reading, i', of a clock which keeps uniform 
time, as v 

a = k x + kg t' + t' (2.4-9) 

where k x and kg are the constant and the rate correction, respectively. 
Equations (2.4-7) to (2.4-9) can serve as a mathematical model to 
compute the right ascension and obliquity p from repeated observa- 
tions of declinations according to the procedure given in Section 2,4.2. 
Consequently, the right ascension will refer to the true celestial 
equator defined by the pole (c'), and the true vernal equinox, i.e. , 
it will be independent of the instantaneous rotation axis (I). In view 
of equation (2.4-7), the solution for the constant Qc' p will include the small 
constant 50 c . ' 

The constant of luni -solar precession is principally derivable from the 
secular variation in right ascension. Detailed procedures for determining 
the constants of precession are found in the literature. 

The constant of nutation, N, which is the coefficient of the 18.6 
year nutation in obliquity is another constant needed in order to be able 
to completely give the true position of (C ' ) in space at any time. 

It is a function of the astronomical constants (A3 - Ai)/A 3, the 
lunar mass, and the obliquity. The constant of nutation determines 
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the common, multiplier (scalar) for all other nutation, terms. Various 
determinations of the constant of nutation are discussed in the literature. 
One of the most complete -determinations is' that by Fedorov [1958]. One 
should note that the observational determination of N includes the small 
term with the same period in Sg^ (equation 2.3-24). This fact has 
to be taken into account when computing the common multiplier of the 
nutations or when one attempts to verify a relation of the type of equa- 
tion (2.2-47). 

2.4.4 Practical Aspects for Adapting a Set of Nutations 

In the usual astronomical work one does not observe the same 
star at upper and lower culmination as discussed earlier; rather, one 
observes stars at either their upper or lower culmination. The latitudes 
resulting from such a scheme are 

§ = 6 ± z (2.4-10) 

in case the upper culmination occurs south or north of the observatory 
zenith. Latitudes computed by this formula always refer to the equa- 
tor of date to which the declinations are counted, i.e., the colatitude 
is the instantaneous angular distance between the observers’ zenith and 
the pole as defined by the set of adopted nutations. According to the 
definition of the celestial pole (C / ) any latitude computed by 

$c' P = 6 c'± z 

varies only due to the progressive Chandler motion. The use of any 
other set of nutations will result in a computable variation of the lati- 
tude. 

Consider the change in declination due to nutation [Mueller, 1969, 
p. 74], 

A6 C / = -{£<•/ ;, sin e c/p COSO! + $./ N sin a (2.4-11) 

This formula is in agreement with our sign convention for l / ) , i.e., 
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positive from west to east, a is the right ascension of the star. 
Equation (2.4-11) can be written as 

A6 c y = Re(i x c /e" lCV ) (2.4-12) 

The difference in latitude based on, for example, the nutations of the 
pole (C^) and the instantaneous rotation axis (I) is 

<£> c ' — j = 6 C / - Si 

_i a 

= He [i (X c / N - Xrw) e ] 

= He [i (x HN + SxS 3 f - Xhn - Sxjf) e ] 

= Re [i (SxS 3 f - 6xr f )e“ ia ] (2.4-13) 

With equation (2.3-20) and (2.3-22) this difference becomes 

$ c / - *x = ^ e" i(AG:j+a) (2.4-14) 

The summation is to be taken over all tidal waves except those 
having sidereal frequency. As was mentioned earlier, the waves of 
sidereal frequency are responsible for the small constant from 59° . 

In case of meridian observations the right ascension O' in (2.4-14) 
can be replaced by the mean sidereal time MST. Thus, from equation 
(2.4-14), the first equation in (2.3-13) and equation (2. 3-15) the latitude 
difference is the negative of the meridional diurnal motion component of 
the rotation axis (I) 

~ $1 — — He (Ujf) = - Uijf (2.4-15) 

The difference is usually referred to as "dynamical variation of latitude ." 
Again, the term with sidereal frequency is excluded in (2.4-15). Equa- 
tion (2.4-14) clearly shows that latitude determinations as derived from 
zenith distance measurements to stars at transit exhibit a nearby di- 
urnal variation. If one observes the same star at consecutive transits, 
the latitude variation shows periods equal to those of the nutations and 
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with, respective amplitudes as given in equation (2.4-14). 

The situation is very similar for longitude observations. The 
usual formula for correcting the right ascensions due to nutation is 

Ag: c > = -d) c ' N cos 0 C / p + sin Q c / p sin<^ - q c / n cosa)tan6 (2.4-16) 

where j/j is again counted positively eastward. The first term, which 
is a correction to the vernal equinox, is called the equation of the 
equinox 

Eq. E — m co ® 9;/ p 

It is equal to the difference between the mean and true right ascension 
of a body on the equator. The other terms, which are a function of 
the position of the star, can be rewritten as 

=Im(-ix c / N e ) tan 6 (2.4-17) 

The basic longitude equation is 

A = MST - GMST 
= AST - GAST 

For reasons of compatability both the local and the Greenwich sidereal 
times are corrected for polar motion. Polar motion, of course, is 
understood to consist only of the motion of the pole (C'). Meridian 
observations give 

A c ' = a - GAST 

= + Eq. E +Im(~iX c , N e _ia )tan 6 - GMST - Eq.E 

— iCt 

= O' H + Im(-i X C , N e ) tan 6 - GMST 

O' H denotes the mean right ascension of date. The longitude computed 
in this manner will vary only due to the progressive Chandler motion. 

Computing the longitude based on the nutations of the instantaneous 
rotation axis (I) gives the, difference * .• 
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—AOL 

A c ' - At = Im[-i(X c ' N - x IN )] e tan 6 

= to e- 1 < A “‘ + “>) tan 6 (2.4-18) 

For transit observations, this difference, according to the first equation in 
(2.3-13) and equation (2.3-15) becomes 

A c ' - Ai = I m (uif) tan 6 = u^ f tan 6 (2.4-19) 

where u' 3 3? is the component of the diurnal motion of the rotation axis (I) 
orthogonal to the local meridian. When verifying the longitude difference 
from observations, one has to account again for the fact that the obliquity, 
as derived from observations, contains the small constant term 5 8°. 

Thus the sidereal terms are excluded in (2.4-18) and (2.4-19). 
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3. MODELS WITH LIQUID CORE 
3.1 Brief Historical Review 

The desire to explain the discrepancies between the observed and pre- 
dicted motion of the earth resulted in a long history of a search for better 
earth models. The first efforts to investigate earth models with liquid core 
date back to Hopkins [1839] . After Chandler’s explanation regarding the basic 
period in polar motion. Hough [1895] and Sludskii [1896] discovered indepen- 
dently the possibility of a nearly diurnal free wobble (NDFW) for a rotating 
container with liquid core. Poincare [1910] published his investigation on the 
precession of the deformable earth, giving an elegant accounting for possible 
movements of the core. His theory is reproduced in [Melchior, 1966]. 

Takeuchi [1950] carried out numerical integration of the equations for a hetero- 
geneous and compressible globe by utilizing different models of the earth’s 
interior, these models being constructed according to seismological results. 
Jeffreys and Vicente [1957 a, b] continued the studies of Sludskii and Poincare 
and proved, based on the work by Takeuchi, that an effect of resonance due to 
the movements in the liquid core (core resonance) appears on the waves whose 
period is close to that of the NDPW. Molodenskii [1961] published an analogous 
study with results similar to those of Jeffreys and Vicente. This theory is 
also reproduced in [Melchior, 1966] . Recently, Shen and Mansinha [1976] pre- 
sented an extension of Molodenskii’s theory. McClure [1976] studies the effect 
of transverse meridional core-flow relative to the mantle on the separation be- 
tween the total angular momentum vector and the rotation vector of the shell 
based on a generalization of Poincare's model. Smith [1977] investigated, 
theoretically and numerically, the free modes based on geophysically plausible 
rotating, slightly elliptical earth models. 
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The above review is necessarily incomplete in view of the large 
amount of work which has been devoted to the subject. The degree 
of difficulties in the mathematical treatment depends on the specific 
assumptions regarding the structure of the core and shell, the core 
boundary , etc. In this section only a discussion of the additional 
characteristic motions is given. For any derivations the reader is 
referred to the literature. 

3.2 Shell and Core Interactions 

The interactions between the core and the shell effect the orienta- 
tion of the shell in space. Determination of the orientation of the 
shell, therefore, can give valuable information about the core motion 
spectrum. 

3.2.1 Free Mode 

The nearly diurnal free wobble (NDFW) is a possible free mode 
which is given by the homogeneous solution of the equations of the 
combined motion of shell and core. It can be categorized as a body (shell)- 
fixed motion of the instantaneous rotation axis (I) having a nearly 
diurnal period. In the astronomical literature it is sometimes refer- 
red to simply as nearly diurnal free "nutation. " Its observational sig- 
nificance has been explained recently. It was Toomre’s [1974] contribu- 
tion to strengthen the point that even this nearly diurnal free "wobble" 
must be accompanied by a free nutation in space whose amplitude and 
period is approximately 460 times larger than that of the NDFW. These 
motions can formally be visualized again by Poinsot's kinematieal repre- 
sentation. The motions occur regardless of whether the shell is taken 
rigid or elastic. The mathematical derivations of this mode can be 
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found in, for example, [Rochester, et al. , 1974] or [McClure, 1976]. 

The frequency of the nearly -diurnal free' wobble, a , is primarily a 
function of the ellipticity of the core as is indicated by the relation 

Ax A3 - a;\ 

which holds for a rigid shell and ideal fluid core. The symbols A, and 
A 3 denote the moments of inertia of the core. A 0 is the equatorial 
moment of inertia of the shell alone, and the other symbols have un- 
changed meaning. The frequency of the accompanying free nutation is 
equal to cl" plus the earth's sidereal motion. This is a property of 
Poinsot's kinematical representation. Thus, the ratio (2.3-29) is 



r _ cl " + Q & 1 

s CL 11 “ 460 


(3.2-1) 


Both the nearly diurnal free wobble and the accompanying 
nutation are retrograde motions as shown in Case 1 of Figure 2,5. 

These notions are super-imposed on those discussed earlier. In par- 
ticular, the shell-fixed motion of the instantaneous rotation axis of the 
shell is 

ic/'t 

u I$ = u x + re (3.2-2) 

The space-fixed motion of the instantaneous rotation axis of the shell 
is a" t 

K ls = X T + 460 re 460 (3.2-3) 

The spatial change, as expressed in equation (3.2-3) is the same for 
■the north celestial’ pole (C / ); thus, a // 

Xr /' = k c ' + 460 r e 460 (3.2-4) 

C" stands for "north celestial pole of the shell . " Probably the most 
meaningful and unique labeling of the two harmonic motions discussed 
here is "retrograde free principal core nutation." The word "core" 
indicates the origin of the motion and the adjective "principal" implies 
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that more motions of this type are possible depending on the assumed 
core structure. Jeffreys and Vicente [1957b], for example, found 
another free but prograde mode because of the density structure in 
their core model. This clearly demonstrates the importance of the 
assumptions about the form and structure of the core -m this type of 
calculation. 

3.2.2, Core Resonance 

Another geophysical phenomenon responsive to the presence of 
the free mode is the tides of the solid earth. Jeffreys and Vicente 
pointed out that nearly diurnal tides whose frequency is close to ot* 
should experience amplifications. This phenomenon is called "core 
resonance. " It is presented mathematically in terms of a factor 
which depends on the tidal frequency and by which the tidal amplitudes 
Aj are multiplied. The factor increases as the tidal frequency ap- 
proaches ol" . Jeffreys and Vicente, as well as Molodenskii, gave es- 
timates for the change in the tidal amplitudes Aj , which can, via equa- 
tion (2.2-43), easily be converted to changes in nutation. An exten- 
sive analysis of these effects based on various models can be found in 
Melchior [1971]. Some of the larger effects are given in Table 3.1. 
The values are to be added to the nutation of the rigid earth. The 
signs are such that the amplitudes of the nutations increase in magni- 
tude. The first column in Table 3.1 denotes the tidal waves in terms 
of their code number, which is explained in Section 2.2. In addition, 
the commonly.. used symbol is given for each tide. The significance 
of these symbols in terms of systematizing the tidal waves is found 
in [Melchior, 1966]. Thus the waves OOi and O x form the semi-monthly 
elliptical nutation. 



Table 3. 1 Corrections for Rigid Earth Nutations 
in Arcsec Due to Core Resonance 


Period (days) 

dj|j sing 

de 

semi-monthly (13.7) 
OOi ; 185.555 

O x 5 145.555 

0.0020 

0.0026 

semi-annual (183) 
(p^; 167.555 

P 1 ; 163.555 

0* 016 

0.020 

annual (365) 
166.554 

S,; 164 556 

-0.0077 

0.0056 


As was mentioned above, it is the relative position of the tidal 
frequency with respect to the frequency of the nearly diurnal wobble, 
cl" , which determines the magnification of the tidal amplitude. cl" 
is not identical to the diurnal frequency of the K x tide but it is closer 
to that of the ^ tide. Since the elliptical nutations consists of two 
tidal waves whose frequencies are symmetric with respect to K x , the 
two waves experience dissymmetric resonance effects. An example is 
the annual nutation which is formed by the tidal waves and S 1 . Both 
tidal amplitudes are of equal magnitude, which results in a zero annual 
nutation in obliquity according to equation (2.2-42). But the liquid 
core model requires a nutation in obliquity solely due to the dissym- 
metric correction. 

The corrections of Table 3.1 have to be added to both the nutations 
of (I) and (C ') . This is so because tidal analysis gives corrections to 
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the tidal amplitudes which enter in both sets of nutations according 
to equations (2.3-23). 

3.3 The North Celestial Pole of the Shell (c") and Its Generalization 

The discussion of the section above demonstrates that for the real 
earth somewhat different terrestrial motions of the instantaneous rota- 
tion axis and a slightly different orientation in space have to be expected 
compared to ideal models, such as the rigid and elastic model or even the 
model with liquid core. Since all geometric observations (directions, 
ranges) take place on the shell one may define the orientation of the 
earth in terms of the orientation of the shell. 

Similar to the definitions of Section 2. 3 . 4, the orientation of the 
shell in space is based on the direction of the pole (C y ), i.e., the 
direction which has neither periodic diurnal motions relative to the shell 
nor to space. A generalization of this definition is obvious. Instead 
of investigating specific models, for an axis which has the desired proper- 
ties of the celestial pole one can simply define the celestial pole with- 
out reference to any model. Thus, 

the orientation of the earth is based on the direction 
of the celestial pole (C), i.e., the axis which has 
. neither periodic diurnal body-fixed nor space-fixed 
motions. The body-fixed motions of (C) are called 
polar motion. 

This definition can be related to the actual measurement procedures. 

One can now say that the observations take place on the surface of the 
earth and that the term "body-fixed" is to be understood as a motion 
relative to a "rigid" surface on which the stations are located. It is 
noted that the effect of- the motion of the crust due to rigid body tides 
has been assumed to be removable by computations. Another type of 
relative station motion, which is very slow and occurs over long periods 
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of time usually referred to simply as ’’crustal motion" or "plate •motion," 
is left out of consideration here. Their implications regarding coor- 
dinate system definition are studied extensively in[-Leick, 197-7].- The 
realization of the direction (C) in the body-fixed frame and in space 
is principally possible through observation of fundamental latitudes 
and declinations (Section 2.4.2). Each term in the set of nutations 
for the pole (C) is empirical and has to be determined from observa- 
tions; whereas, the relative magnitude of the nutation coefficients for 
a specific model result from the motion theory. Fortunately, the rigid 
or the elastic model approximates the actual motion of the earth very 
well. It is, therefore, not necessary to attempt to observe the whole 
spectrum of nutations; rather, one can limit the investigations to specif- 
ic frequencies whose amplitudes are most likely to deviate from the model 
values. The frequencies of Table 3.1 are certainly among such "danger- 
ous" frequencies. 

Formally, one can write the space-fixed position of the celestial 
pole (C) as 

He = He' + Ahfcn + Ah c r + A Xu (3.3-1) 

where x c ' (equation 2.3-25) can be considered the first order approxi- 
mation of the pole (C); whereas, the remaining empirical terms are of 
second order. The term Ah Fcn denotes the contribution of the free core 
nutation. This term is entirely empirical. Neither its magnitude nor 
its phase is known since it is a result of the homogeneous solution. The 
term Ax CR denotes the resonance effect on the nutation due to core 
motions. In this case, the frequencies and estimates of their amplitudes 
are obtained from theory (Table 3.1), but depend strongly on the assump- 
tions of the core structure. Therefore, these terms are essentially 
also of empirical nature. The last term, Axu> stands for all effects 
which are not yet accounted for. 
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In the usual latitude observation work where a star is observed 


only at one culmination, as represented by equation (2.4-10), 

= 6 ± z 


the computed colatitude refers to the axis defined by the adopted set of nuta- 
tions, This axis has no physical meaning. It only serves to transform the 
actual, observable into a latitude -like quantity. Thus, the additional shell- 
fixed motion of the instantaneous rotation axis (I) which is due to the nearly 
diurnal free wobble (equation 3.2-2) is irrelevant for this observational proce- 
dure. Therefore, regarding the determination of the nearly diurnal free wobble, 
the only thing one can do is to determine the associated free nutation from the 
analysis of the observed zenith distance and conclude then via equation (3. 2-1) 
the actual size of the nearly diurnal free wobble . Using equations (2 . 4-10), 
(2.4-12), and (3.3-1), the computed latitude variation will be 


A# = $ c - # c / 

— Re i(Axpc N + Axcr + A Xu) e ^ 

(3.3-2) 

= ASfcn + A6cr + A6u 
Ulc diurnal 

In principle, any adopted set of nutations can serve in this type of 
analysis. Instead of taking equation (2.3-25) for the pole (C 7 ), one could 
just as well select the nutations of the rotation axis (I) of (2.3-20) or the 
angular momentum axis (H) of (2.2-43). In these cases, additional 
computable terms would appear on the right-hand side of equation (3.3-2). 

The corrections in declination in equation (3.3-2) are identical to the 
negative of the component of the diurnal motions of the adopted pole (C') around 
the pole (C) along the local meridian. It is emphasized that this diurnal motion 
strictly refers to the pole (C') as defined by the adopted set of nutations . Ithas 
nothing to do with the diurnal body-fixed motion of the actual instantaneous ro- 
tation axis (I) or angular momentum axis (H) . These latter motions are not 
observable! 



It is realized that the latitude variation of equation (3.3-2) is very slow 
if zenith distances are observed at culmination of the same star. It is, there- 
fore, permissible to approximate the diurnal motion by a circle of constant 
radius over an interval of, say, one day, five days, etc. A <3? has then the 
same magnitude for each station's observation. This difference is part of the 
Kimura[1902] term which was introduced in order to represent the "non- 
polar variation in latitude. " The Kimura term also absorbs other constant 
effects which are specific to astronomical observation, such as errors in 
proper motion, aberration, etc. In any case, it is important to note that the 
introduction of Kimura’s term makes it, at least conceptually, possible to 
determine the celestial pole (C) even from observations at one culmination 
only. One only needs the zenith distance observations from several stations 
of the same star at culmination in order to determine polar motion and the 
constant non -polar term in the least squares sense. Let ‘J’u be the adopted 
latitude, then the equations are in the usual notation: 

$u = ®c-«iccos A - u 2 c sin A 

- q/ — m c cos A — u^ c sin A + A6 FCN + A6 C r + A 5u 

This expression can be rewritten as 

- ®c' = cos A - Uac sin A + Z 

where 

Z Ulc/ diurnal 

is the Kimura term. The subscript c' denotes the adopted set of nutations . 

Since in all computations, the rigid earth nutations of (I) have been 
used, the analysis of Kimura’ s term should give information on core reso- 
nance, etc. In fact, the International Latitude Service introduced in 1955 a 
special observation program, which is called the three-group observations, 
in order to make an effective analysis of the Z-term possible. 
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3.4 Observational Evidence - IAU (1977) Proposed Nutation Series 

Reports claiming to have observed the nearly diurnal free wobble, 
NDFW, based on above or below pole observations, were given at 
several occasions. A summary is given in Yatskiv [1972]. The earliest 
observation, related to this second free mode, seems to have been re- 
ported by Popov [1963]. He found an amplitude of 0". 016 from a long 
series of latitude observations of two stars culminating approximately 10 
hours apart. Until the clarification made by Toomxe [1974] the astron- 
omers, apparently, being unaware of the associated large free nutation, 
interpreted their observations indeed as those of the NDFW. Exam- 
ining the method of analysis of Popov along the lines discussed in pre- 
vious sections, it becomes immediately clear that the x'eported variations 
actually represent the change in nutation. The latest report on the NDFW 
is [Yatskiv et al. , 1975]. They not only confirm the presence of the 
retrograde mode, but also find that the existence of the prograde mode 
is quite possible on the basis of the available data. 

Some of the effects of core resonance have also been confirmed 
by observations. A correction of 0". 02 for the semi-annual nutation 
was obtained from analysis of Kimura’s term [Wako, 1970]. The 
presently adopted constant of nutation for the epoch 1900.0 is N - 9". 21. 
Using the current best estimates of the mass of the moon and the con- 
stant of precession, the relation (2.2-48) yields a somewhat larger value 
of N = 9 n ,22. But from observations one arrives constantly at a smaller 
value N = 9 n . 20 [Fedorov, 1958]. Applying the correction due to core- 
resonance to the computed value reduces the latter to the observed value. 
This is certainly one of the nicest features of the liquid core models. 

It explains one of the longest known incohereneies in the system of 
fundamental constraints., 
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An astronomical determination of the fortnightly term was reported 
by Guinot [1970] for longitude observations and by Morgan [1952] for 
latitude observation. Morgan’s analysis was based on P.Z.T. observa- 
tions at Washington. While using the adopted set of nutations of the 
rotation axis (I) for the rigid earth, he obtained a fortnightly correction 
equal to that of the Oppolzer term, just as rigid model theory predicts. 
More complete analyses were carried out by Fedorov [1958] in terms of 
the fortnightly diurnal motion term, and recently by McCarthy [1976], who 
computed the corrections to the fortnightly nutation directly. Both ob- 
tained values which are in agreement with those predicted by the co- 
resonance model (Table 3.1). They used the nutations of the angular 
momentum (H) as computed from Poisson's equation, as their reference 
pole. One may add that McCarthy used (H) as an approximation to the 
rotation axis of the rigid model and corrected the observations for the 
Oppolzer terms (equation 2.4-15); whereas, Fedorov specifically intended 

I 

to use (H) because its position in space is relatively independent of 
mass redistribution. But from the merely analysis point of. view, this 
distinction is unimportant since the adopted set of nutations (Poisson 
solution) serves only as an intermediary reference standard. McCarthy’s 
corrections should be interpreted as corrections to the nutations of {C r ) 
in order to get the nutation of (C). Fedorov's corrections need to be 
added to the nutation of (H) in order to get the nutation of (C). This 
is so because astronomical observations give the motion of the pole, as 
defined by the adopted set of nutations, with respect to (C). The 
fortnightly correction term of Table 3 . 1 can be readily converted to 
a corresponding diurnal body-fixed motion term. The derivation of its 
radius follows from equation (2.2-46) 

H 3 s “ a 3S cos (2 s) + i b 2s sin (2s) (3.4-1) 

where 2 s denote's the fortnightly term, and from equation (2.3-19), 
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(3.4-2) 


Ugs — i A h 2e 



Neglecting the subscripts and combining equations (3.4-1) and (3. 4-2), one gets 


Ua £ 


a + K Q „ "K 

— - — sin (<p - 2s) + — - — sin ((f) + 2s) 


- l 


a ^ cos ((f) - 2 s) + i cos (<p + 2 s ) 


2 ~~ - 2 

(3.4-3) 

Substituting the coefficients from Table 3.1 and neglecting their differences 
gives 

6u 2s = -0.0023 sin (<p - 2s) - i 0.0023 cos(<p-2s) (3.4-4) 


Fedorov found a radius for the fortnightly term of 0"009. As far as 
the model is concerned, this radius formally corresponds to the rigid earth 
nutation plus the core resonance effects, i.e., as predicted by Jeffreys and 
Vicente [1957, a, bj. Equations (2.3-17) and (3.4-4) give 

|uas| = 0.0052 + 0.0023 = 0’.'0085 

The perfectly elastic model, however, predicts a fortnightly diurnal term of 
2 

approximately ~ 0'.'0062, as is seen from the second equation in (2.3-13). 

We may therefore conclude that, at least for the fortnightly nutation, the 
perfectly elastic model does not seem to conform with observation. Since 
the BIH [1975] has currently adopted the nutation of the pole (C 7 ) of the elastic 
model (equation (2. 3- 25)), their procedure may require revision at some 
future time. 

It appears, based on the results of Table 3.1 and the method implied, 
that the actual diurnal radius of the motion of (E) around the pole (C) does not 
change significantly because of core resonance effects. This statement is 
possible because it can be shown that the same corrections to the tidal ampli- 
tudes cause nearly identical corrections for both (I) and (C 7 ) (see equations 
2.3-23). However, this fact cannot be proven from astronomical observations 
since (I) is not observable. 
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The IAU-Symposium No. 78 (1977) recommended that the following 
set of coefficients be substituted for the corresponding coefficients in 
Woolard’s series for the nutations in-order to provide 'a more accurate 
representation of the forced nutation of the axis of rotation of the 
earth due to the luni -solar perturbing forces: 

Table 3.2 IAU (1977) Proposed Nutation Coefficients 


Period (days) 

0 sin 0 

0 

6798 

-6.843 

9.206 

3399 

0.083 

-0.091 

365 

0.058 

0.006 

183 

-0.520 

0.569 

122 

-0.020 

0.022 

27.6 

0.028 

0.000 

13.7 

-0.083 

0.091 


The sign convention in Table 3.2 is the same as in the AENA Supplement 
[1961, p. 44] . The coefficients are obtained by adding the tidal cor- 
rections to the nutation of the instantaneous rotation axis (I) of the 
rigid body. The thus modified nutation series does not describe the 
position of the celestial pole (C). The theoretical corrections for the 
dynamical variations are still needed! 
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4. SUMMARY AND RECOMMENDATIONS 

The motion characteristics of the rigid, the elastic earth model 
and the model with liquid core have been reviewed. An attempt has been 
made to strictly distinguish between motions due to external forces and 
the free motions. The dynamical theories of the rigid and elastic earth 
model involve various axes, such as the axis of rotation (I), the angular 
momentum axis (H), the axis of figure (P) and an axis called the celes- 
tial pole (C / ). Special efforts have been made to investigate the sig- 
nificance of these axes regarding observability. It has been found that 
even for these ideal models, where the stations do not change due to 
crustal motions, only the celestial pole (C r ) is observable through 
fundamental astronomical observations. "Observable" is to be under- 
stood in the sense that the direction between the observatory and the 
celestial pole (C / )can be measured without using any hypotheses or models. 
The pole {C^ moves with respect to the body only because of the pro- 
gressive Chandler motion. It is understood that fundamental astro- 
nomical observations, where the same star is observed at both cul- 
minations, give only the mean positions for the twelve hour time span. 

The diurnal motions of the instantaneous rotation axis (I), and the 
angular momentum (H) have been investigated. The periodic diurnal polar 
motions of (I) are strictly related to the spatial nutations of (I). The 
correspondence between those motions can be demonstrated most easily 
in terms of two cones rolling on each other, having characteristic 
vertex angles and speeds. It has been shown that the instantaneous 
rotation axis (I) is needed only at the initial step when formulating the 
motions. From the descriptive point of view, one can associate with 
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a rotating body which is subject to external forces two instantaneous 
rotation axes. One of the two axes is (I), which must have periodic 
diurnal body-fixed motions. The other axis is (C*); it is effected 
only by the progressive motion (free motion). Since the Chandler 
(Euler) motion, at least conceptually, represents the initial conditions, 
it can assume any magnitude. (Please remember that we are talking 
about ideal bodies.) If its magnitude is zero, then the celestial pole 
(C') does not move with respect to the model surface; whereas, (I)- still 
has the diurnal motions around (C / ). This type of explanation is 
strictly valid since the forced and free motions are independent, as 
they represent two independent solution components to the differential 
equation of motion. As for adopting a set of nutations, clearly, pref- 
erence has to be given to the celestial pole C*, not only because it is 
observable (in case of rigid and elastic bodies) but also because the 
positional elements which refer to it have no diurnal periodic variations 
and do not require a correction for the so-called dynamical variations. 

The differences and commonalities between the rigid and perfectly 
elastic model have been discussed at length, the main characteristic 
being that the angular momentum axis remains virtually unchanged in 
space whereas the direction of the (C') changes somewhat. But it 
was pointed out that the observational evidence does not entirely confirm 
the predicted changes. 

The liquid core model is the most general model considered in 
this study. It introduces a new spectrum of motions, such as the nearly 
diurnal free wobble (NDFW) and its associated change in nutation, both 
resulting from the free solution. In addition, certain frequencies change 
their amplitudes due to core resonance. The observations fit such a 
model quite well. Similarly, as in the case of the rigid, and elastic 
bodies, the celestial pole of the shell ( C ") was defined as that pole 
having neither shell-fixed nor space-fixed periodic diurnal motions. Thus, 
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in the case of the liquid core model, the celestial pole of the shell (C*) 
is observable. 

Above it has been said that the celestial pole (C*) or (C /7 ) is 
observable depending on the mathematical model under consideration. 

The real earth, of course, does not behave exactly as these models in- 
dicate although the liquid core model fits the observations better than 
the other two models. Yet, the concept of the celestial pole, or even 
better that of the celestial pole of the shell can be extended so as to 
denote that pole which is observable in actuality. This pole can then 
simply be called the "Celestial Pole (C)." It has the property of having 
no body-fixed and n© space-fixed diurnal periodic motions. The naming of 
this observable pole correctly does not give any hint as to the best 
fitting theoretical model since, anyhow, each model is only an approxi- 
mation. The point of view is taken that the nutations and, of course, 
the polar motions of the celestial pole (C) can only be determined from 
observations. The possible adoption of the IAU - 1977 set of nutations 
is a step in this direction. Unfortunately, this set of nutations does not 
include the theoretical Oppolzer terms; therefore, the dynamical cor- 
rections are still needed. 

It is recommended that in the future the terms "polar motion" and 
"nutation" be only associated with the body-fixed and space-fixed 
motions of the celestial pole (C), respectively. If the motion of any 
other axis is meant, the name of this axis should be given explicitly. 

The celestial pole (C) gives the natural reference direction not only 
for astronomical observations but also for laser ranging to lunar reflec- 
tors. More details on this subject are given in[Leick, 1978], 

Finally, it is sometimes suggested that the nutations be given for 
the angular momentum (H). But it is clear that such a procedure violates 
the concept of observability. It may be true that the direction of (H) in 
space is the same for the rigid, elastic, or any other reasonable earth 
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model. But this property is not of much interest to the astronomer 
or geodesist who tries to determine the orientation of the earth. It 
is conceptually simpler to refer to an axis which is observable. 

Both representations are actually equivalent, i.e. , determining the nuta- 
tions of the celestial pole (C) or using the adopted (rigid model) nuta- 
tions of (H) but then determining the periodic diurnal body -fixed motions 
of the angular momentum (H). 
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